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Measurements  of  optical  and  paramagnetic  resonance  spectra  made  on 

yttrium  aluminum  garnet  doped  with  0.1  percent  chromium  have  been  suc- 

4 
cessfully  fitted  to  existing  theory.   The  locations  of  the  FT  and 
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Fr  ,   Pr   levels  were  determined  to  be  16,700,  23,200,  and  35,700  cm   , 

respectively,  by  optical  absorption  measurements.   From  these  results, 

the  crystal  field  parameter  Dq  was  found  to  be  1670  cm   .   The  locations 

of  the  two  T  levels  were  fitted  within  3.9  percent  to  the  theory  of 

4        4 
Finkelstein  and  Van  Vleck  which  predicts  that  the  FT     and  PT   states 

are  coupled  by  the  octahedral  crystal  field  perturbation.   The  results 
of  paramagnetic  resonance  measurements  give  the  g- factor  =  1.980  and  the 
zero  magnetic  field  ground  level  splitting  =  0.509  cm  .   These  results 
were  used  to  determine  the  splitting  of  the  excited  T     level  due  to  the 
trigonal  distortion,  which  is  found  to  be  2140  cm  .   The  locations  of 
the  two  component  energy  levels  of  the  V     level  are  then  used  to  calcu- 
late the  contributions  to  the  crystal  field  potential  expansion  caused 


v 


by  the  trigonal  distortion  by  a  method  which  treats  the  distortion  as  a 
small  perturbation  of  the  T      level.   This  perturbation  calculation  was 
carried  out  to  second  order  to  determine  the  series  expansion  coeffi- 
cients of  this  distortion.   Only  the  coefficient  of  the  second  order 
spherical  harmonic  y   is  found  to  contribute  to  the  deviation  from  a 
pure  octahedral  cubic  crystal  field.   The  corrected  potential  expansion 
was  then  used  to  determine  the  splittings  of  the  two  T      levels . 
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SECTION  I 
INTRODUCTION 

I. A.   Brief  Historical  Sketch 
Paramagnetic  resonance  may  be  observed  in  almost  any  substance  which 
has  a  net  electron  spin  different  from  zero.   Examples  are:   organic  free 
radicals,  some  free  atoms,  and  ions  with  a  partly  filled  electron  shell. 
In  many  instances,  paramagnetic  resonance  (PMR)  is  observed  in  substances 
contained  in  host  crystals.   In  materials  of  this  type,  the  PMR  absorption 
is  actually  caused  by  the  dopant  ion,  but  the  structure  of  the  host  crys- 
tal contributes  greatly  to  the  PMR  behavior.   For  this  reason,  PMR  has 
proven  to  be  a  valuable  method  for  the  study  of  the  crystalline  field,  es- 
pecially when  complemented  by  optical  absorption  measurements.   Ions  of  the 
iron  and  rare  earth  groups  have  been  studied  extensively  by  this  technique. 

The  first  magnetic  resonance  measurements  were  performed  by  E.  Zavoisky 
in  Russia  in  1945.   He  used  a  radio  frequency  source  with  a  wavelength  of 
25  meters  to  excite  resonance  transitions  in  several  substances  at  very  low 
magnetic  fields,  at  which  clearly  defined  resonances  were  scarcely  discern- 
ible.  The  availability  of  sophisticated  microwave  sources  subsequent  to 
World  War  II  had  a  great  impact  on  this  area  of  experimentation,  and  Z avoisky 
and  other  workers  used  these  sources  to  expand  PMR  measurements  into  the 
microwave  frequency  region.   Much  of  the  early  work  was  done  by  the  groups 

at  Oxford  who  have  made  most  of  the  contributions  to  both  the  PMR  theory  and 

1,2 

measurement  techniques.   Their  many  publications  '   provide  a  valuable  store 

of  information  for  the  study  of  PMR  and  crystal  field  theory.   During  the 
past  few  years,  PMR  has  been  observed  in  the  optically  excited  states  of 

chromium  and  manganese  by  a  double  resonance  technique  developed  by  Geschwind 

3 
and  others . 


In  addition  to  the  spin-orbit  interaction  perturbation  which  acts 
upon  free  ions,  dopant  ions  in  crystalline  hosts  are  also  acted  upon  by 
the  crystal  field  perturbation.   For  iron  group  ions,  this  perturbation 
is  stronger  than  the  spin-orbit  interaction,  and  the  various  energy  levels 
resulting  from  the  crystal  field  are  usually  separated  by  energies  which 
may  be  determined  by  optical  absorption  measurements.   For  this  reason, 
these  measurements  may  be  used  as  the  basis  for  calculating  the  coefficients 
in  the  expansion  of  the  crystal  field  potential .   PMR  measurements  comple- 
ment optical  measurements  in  the  calculation  of  the  spin-orbit  interaction 
coefficient  and  smaller  crystal  field  distortions,  such  as  departures  from 
high  symmetry,  if  present.   For  rare  earth  ions  in  crystalline  hosts,  the 
spin-orbit  interaction  is  stronger  than  the  crystal  field  perturbation  be- 
cause the  incomplete  inner  shell  of  electrons  is  partially  shielded  by  an 
outer  shell.   In  this  case,  the  spin-orbit  states  are  separated  by  energies 
measurable  by  optical  absorption  techniques,  and  results  of  optical  measure- 
ments may  be  used  to  calculate  the  spin-orbit  coupling  parameter .   The  re- 
sults of  PMR  measurements  may  then  be  used  to  calculate  the  crystal  field 

4  5 
parameters,  as  was  done  by  Elliott  and  Stevens  '   for  several  rare  earth 

ethyl  sulfates.   Complementary  optical  and  PMR  absorption  measurements 

therefore  provide  a  means  of  characterizing  the  effects  of  crystal  field 

and  spin-orbit  perturbations  on  the  energy  levels  of  both  iron  group  and 

rare  earth  group  ions  in  host  crystals. 


I.B.   General  Crystal  Field  Theory 

The  effects  of  the  crystal  field  on  the  energy  levels  of  dopant  ions 
in  host  crystals  are  evaluated  by  expanding  the  crystal  field  potential 
V(r,0,4>)  in  a  series  of  solid  harmonics 

V(r,e,<t>)  =  E  A£m  r£  y^  (0,<t>), 
£  ,m 

where  the  A   are  the  expansion  coefficients,  y   are  spherical  harmonics, 
and  r,0,(()  are  the  coordinates  at  which  the  potential  is  measured.  The 
effects  of  the  crystal  field  on  the  energy  levels  of  the  dopant  ion  are 
determined  by  calculating  the  matrix  elements  of  the  above  potential  using 
the  free  ion  states  which  are  eigenfunctions  of  the  wave  equation  including 
the  coulomb  interaction  between  electrons,  but  not  including  the  spin-orbit 
interaction.   This  matrix  is  then  diagonalized  to  give  the  eigenvalues  of 
the  total  Hamiltonian  which  are  then  substituted  into  the  original  crystal 
field  matrix  equation  to  obtain  eigenfunctions  diagonal  with  respect  to 
this  Hamiltonian.   It  may  be  convenient  to  perform  the  calculation  described 
above  for  a  high  symmetry  crystal  field,  and  to  treat  slight  departures  from 
this  high  symmetry  as  perturbations  of  the  states  resulting  from  this  cal- 
culation.  This  approach  is  useful  in  the  study  of  the  trigonal  distortion 
present  in  several  iron  group  ions  in  various  crystalline  hosts.   The  spin- 
orbit  interaction  and  magnetic  field  perturbation  are  then  considered  as 
perturbations  of  the  crystal  field  states  determined  in  the  calculations 
described  above.   The  true  crystal  field  energy  levels  are  determined  by 
optical  absorption  measurements ,  and  the  measured  values  of  energy  are 
fitted  to  the  theoretical  values  by  varying  the  expansion  coefficients  A   . 


The  object  of  the  crystal  field  calculations  is  therefore  to  obtain  values 
of  the  expansion  coefficients  which  give  the  best  fit  to  measured  data  when 
substituted  into  the  crystal  field  potential  expansion. 

I.C.   General  PMR  Theory 

Consider  an  ion  of  mass  m  possessing  a  resultant  angular  momentum  J, 
placed  in  a  constant  external  magnetic  field  "?f.   The  angular  momentum  vector 
of  the  ion  then  precesses  about  the  axis  of  H  with  the  angular  frequency 

~lu  =  g  (~-)"lf,  (1) 

2mc 

where  g  is  the  spectroscopic  splitting  factor,  e  is  electronic  charge,  and 
c  is  the  speed  of  light.   Now  assume  that  a  circularly  polarized  radio  fre- 
quency magnetic  field  is  applied  to  the  ion  with  such  a  sense  and  frequency 
that  it  rotates  about  H  in  synchronism  with  the  angular  momentum  vector. 
This  vector  will  then  ultimately  reverse  its  projection  on  H  with  a  resultant 
exchange  of  energy  with  the  radio  frequency  radiation  field.   This  phenomenon 
is  called  paramagnetic  resonance. 

The  magnetic  moment  u  of  the  ion  is 

tr  =  "#  ,  (2) 

ft 
in  which  3  is  the  Bohr  magnetron  and  ft  is  Planck's  constant.   The  energy  of 
the  ion  in  the  magnetic  field  is  then 

E  =  -tr-H.  (3) 

Taking  the  magnetic  field  direction  along  the  Z  axis  and  substituting 
(2)  for  "y  gives  for  the  expectation  value  *\E>  of  the  energy 

<E>  =  h1  Hz  <Jz>  =  90Hm'  (4) 


where  m  is  the  magnetic  quantum  number.   In  the  magnetic  field,  the  energies 
of  the  ion  corresponding  to  the  various  spatial  orientations  of  J  are  then 
g$Hm,  and  the  selection  rules  for  the  allowed  transitions,  which  are  magnetic 
dipole  in  character,  are  Am  =  ±1.   Transitions  between  the  various  levels  for 
which  this  selection  rule  holds  may  be  induced  by  quanta  of  energy 

fioo  =  gBH,  (5) 

which  is  seen  to  be  identical  with  Equation  (1) . 

Consider  now  a  system  of  the  ions  discussed  above  situated  in  a  crystal 
lattice  and  in  thermal  equilibrium  with  their  surroundings.   These  ions  may 
be  considered  distinguishable,  and  the  ionic  energy  levels  in  the  magnetic 
field  will  be  populated  according  to  Boltzmann  statistics.   At  high  tempera- 
tures, the  various  states  will  be  nearly  equally  populated,  and  there  will 
be  little  net  energy  absorption  at  resonance,  because  transitions  from  a 
lower  to  a  higher  state  and  transitions  from  a  higher  to  a  lower  state  are 
almost  equally  probable.   The  strength  of  PMR  absorptions  is  therefore  en- 
hanced at  lower  temperatures  because  the  population  difference  is  greater, 
resulting  in  a  greater  energy  absorption  by  the  ions.   This  resonance  ab- 
sorption can  be  detected  by  measuring  the  energy  lost  from  the  radio  fre- 
quency field,  which  causes  a  damping  of  the  tuned  circuit  in  which  the  para- 
magnetic sample  is  placed. 

I.D.   Outline  of  Procedure 


Both  optical  and  paramagnetic  resonance  spectra  of  the  Cr    ion  in 
various  crystalline  hosts  have  been  extensively  studied.   Much  of  the  work 

has  been  done  on  the  various  chrome  alums  by  Finkelstein  and  Van  Vleck, 
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Kleiner,   Bleaney,   and  Davis  and  Strandberg.    Other  host  materials  have 


been  studied  by  Low   who  treated  MgO,  Zaripov  and  Shamonin,   who  examined 

12 
ruby,  and  Carson  and  White,  "  who  studied  yttrium  aluminum  garnet  and  yttrium 

gallium  garnet.   Much  of  the  interest  in  Cr  +  in  recent  years  has  been  oc- 
casioned by  the  potential  of  this  ion  as  a  laser  material,  and  its  application 
as  an  active  ion  in  the  ruby  laser. 

The  crystalline  host  yttrium  aluminum  garnet  (YAG)  has  found  broad  ap- 
plication as  a  host  material  for  laser  ions,  especially  the  rare  earth  ions 

.^3+     3+        3+ 

Nd   ,  Er   ,  and  Ho   .   Since  YAG  is  a  hard  crystal  with  high  thermal  con- 
ductivity, and  is  available  in  large  crystals  of  excellent  optical  quality, 

13 
it  is  well  suited  to  be  a  laser  host  material.    Although  it  has  been  deter- 

14 
mined  by  Burns  et  al.    that  chromium  doped  YAG  is  an  inferior  laser  material 

to  chromium  doped  aluminum  oxide,  e.g.  ruby,  the  material  is  of  interest  be- 
cause of  its  similarity  to  ruby  and  because  the  basis  for  the  crystal  field 
and  paramagnetic  resonance  (PMR)  calculations  may  be  used  for  other  materials. 

Trivalent  chromium  has  three  d  electrons  and  its  ground  state  is   ?,,«. 

12 
In  the  YAG  crystal,  it  occupies  a  site  of  octahedral  symmetry   being  sur- 
rounded by  an  octahedron  of  six  oxygen  ions .   Superimposed  on  this  octahedral 
field  is  a  small  trigonal  distortion  caused  by  a  slight  displacement  of  the 
Cr    ion  along  the  local  and  [ill]  crystal  direction.   This  distortion  ac- 
counts for  the  interesting  optical  and  magnetic  properties  of  the  material. 
Figure  1  shows  the  evolution  of  the  energy  levels  in  YAG:Cr  under  the  in- 
fluence of  the  perturbations  to  be  considered  in  this  work.   The  perturbations 

are  shown  in  the  order  in  which  they  will  be  treated.   The  octahedral  crystal 

4 
field  splits  the  Fo/;?  ground  state  of  the  free  chromium  ion  into  a  singlet 

ground  state  T      and  two  higher  lying  triplets  [T    ,    T    ).   Each  of  these  levels 
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is  fourfold  degenerate,  including  Kramers  degeneracy,  so  that  the  total 

4 
number  of  states  in  the  F  ground  manifold  is  4-7  =  28.   Each  of  the  higher 

lying  triplet  states  is  further  split  into  two  crystal  field  sub-levels  by 

the  action  of  the  trigonal  distortion,  in  agreement  with  the  predictions  of 

group  theory.     Second  order  spin-orbit  coupling  interaction  between  the 

ground  (I^)  level  and  the  first  excited  triplet  (V    )  level,  together  with 

the  trigonal  distortion,  work  to  split  the  ground  level  into  two  Kramers 

doublets,  characterized  by  m  =  +3/2  and  m  =  ±1/2.   Upon  application  of 

o  s 

a  magnetic  field  H,  paramagnetic  resonance  transitions  may  be  observed  be- 
tween these  various  doublet  levels,  providing  a  means  of  determining  the 
splitting  of  the  ground  state  doublet  levels  and  the  spin-orbit  coupling 
parameter. 

The  effect  of  the  trigonal  distortion  on  the  T      crystal  field  level 
may  be  determined  by  treating  this  distortion  as  a  small  perturbation  on 
the  larger  octahedral  cubic  field.   The  perturbation  splits  the  triply 
degenerate  V      level  into  a  doubly  degenerate  level  and  a  non-degenerate 
level.   The  magnitude  of  this  splitting  is  determined  by  a  calculation 
based  on  the  PMR  measurements.   The  knowledge  of  this  splitting,  together 
with  the  center  of  gravity  condition  for  energy  levels,  enables  one  to  de- 
termine the  location  of  the  T      sublevels ,  from  which  the  constants  in  the 

b 

crystal  field  potential  expansion  contributing  to  the  trigonal  distortion 
may  be  determined.   Knowledge  of  these  constants  provides  a  means  of  deter- 
mining the  effect  of  the  trigonal  distortion  on  all  of  the  crystal  field 
sublevels. 


The  upper  triplet  crystal  field  level  (T  )  interacts  with  a  higher 

4 
lying   P  level  through  the  crystal  field  interaction,  causing  the  T   level 

4 
to  be  depressed  and  the   P  level  to  be  raised  in  energy,  with  appropriate 

mixing  of  the  two  states.   The  theory  of  this  interaction  has  been  worked 

out  in  detail  by  Finkelstein  and  Van  Vleck,   who  predict  the  location  of 

all  of  the  crystal  field  levels  in  chrome  alum  in  terms  of  the  crystal 

field  parameter  Dq.   Optical  absorption  measurements  made  on  chromium  doped 

YAG,  to  be  described  later,  agree  well  with  this  theory. 


SECTION  II 
THEORY  OF  THE  CRYSTALLINE  FIELD  IN  YAG  DOPED  WITH  CHROMIUM 

II. A.   Expansion  of  Crystalline  Field  in  Solid  Harmonics 

The  crystal  structure  of  YAG,  Y  Al  (A10  )  ,  is  body-centered  cubic 
and  the  unit  cell  contains  eight  molecular  formula  units  or  160  ions . 
Each  of  the  metal  ions  is  surrounded  by  a  polyhedron  of  oxygens,  which 
is  a  distorted  cube  for  the  yttrium  ions,  and  is  either  a  tetrahedron 
or  an  octahedron  for  the  aluminum  ions.  All  experimental  evidence  ob- 
tained thus  far  indicates  that  the  chromium  dopant  ions  replace  aluminum 
only  in  octahedral  sites.   The  chromium  ion  may  therefore  be  viewed  as 
being  at  the  origin  of  a  cartesian  coordinate  system,  with  six  oxygen 
ions  placed  at  equal  distances  from  the  origin  on  the  ±X,  ±Y,  and  ±Z 
axes.   Superimposed  on  the  octahedral  field,  there  is  a  small  trigonal 
distortion  caused  by  a  slight  displacement  of  the  chromium  ion  along  the 
local  and  crystalline  \lll\    direction.   It  is  the  objective  of  the  study 
of  the  crystal  field  to  determine  the  effect  of  the  octahedral  field  and 
trigonal  distortion  on  the  energy  levels  of  the  chromium  dopant  ion. 

Assuming  that  the  crystal  field  potential  V.  satisfies  Laplace's 

equation,  it  can  be  expanded  in  spherical  harmonics  Y   (6,<t>)  as  follows: 

V  =  Z  V . ,  with 
J   D 

V  =  E  A™   /  y£m  (6,*),  (6) 

3   £,m 

where  the  A   are  expansion  coefficients  and  r  is  the  radius  vector  to 

£   m 
the  point  at  which  the  potential  is  evaluated.   The  product  r  y    (6  ,cj>) 

is  a  solid  harmonic.   The  V.  notation  is  used  to  indicate  that  the  potential 

operates  only  on  the  wave  function  of  the  j    electron,  even  though  the 

potential  is  identical  for  all  electrons  outside  closed  shells. 
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In  evaluating  the  effects  of  the  potential  (6)  on  the  Cr   ion,  the 
fact  that  the  crystal  field  energy  level  splittings  within  a  given  multiplet 
may  be  comparable  in  magnitude  to  the  splittings  between  multiplets  caused 
by  the  coulomb  interaction  between  electrons  must  be  considered.   In  order 
to  be  strictly  correct,  the  coulomb  interaction  and  crystal  field  perturba- 
tions must  be  simultaneously  applied  to  hydrogenic  wave  functions,  and  the 
resulting  matrix  must  be  diagonalized  numerically.   This  calculation  has 
been  made  by  Sugano  and  Tanabe,   who  treated  several  transition-metal  ions 
by  this  method.   Finkelstein  and  Van  Vleck  have  determined  the  energy 
levels  of  potassium  chrome  alum  by  a  different  method.   They  first  determined 
the  matrix  elements  of  the  crystal  field  potential  within  each  of  the  multi- 
plets using  the  free  ion  Russell-Saunders  states  as  basis  functions.   The 
resulting  matrices  were  then  diagonalized  individually  for  each  multiplet 
level,  and  a  new  set  of  basis  functions,  simultaneously  diagonal  within  a 
given  multiplet  with  respect  to  the  free  ion  Hamiltonian  and  the  crystal 
field  perturbation,  were  calculated.   This  new  set  of  basis  functions  is 
called  the  LT  representation.   The  matrix  elements  of  the  crystal  field 
potential  were  then  evaluated  between  the  various  multiplet  levels  using 
this  Lr  representation  as  a  basis.   To  the  matrix  of  the  crystalline  poten- 
tial calculated  in  this  way,  the  diagonal  matrix  of  the  Russell -Saunders 
energy  was  added.   The  resulting  matrix  was  then  diagonalized  to  give 

energy  levels  of  the  chromium  ion  in  the  crystal  field.   This  procedure 

4      4 
was  followed  in  this  work  to  determine  the  energy  levels  of  the  F  and  P 

multiplets  of  chromium  in  the  YAG  crystal,  and  will  be  discussed  in  some 

detail  in  the  following  sections. 
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The  number  of  terms  in  the  expansion  (6)  which  contribute  to  the 
potential  can  be  considerably  reduced.   The  term  for  1=0,   which  is  a 
constant  term,  can  be  dropped  because  it  has  no  effect  on  the  relative 
energies  of  the  crystal  field  states.   Furthermore,  all  terms  for  which 
I  >   4  will  have  zero  matrix  elements  and  can  be  omitted,  because  in 
evaluating  integrals  of  the  form  /u*Uvdt ,  where  u  and  v  denote  d-electron 
wave  functions  and  U  is  a  spherical  harmonic,  the  density  u*v  does  not 
contain  harmonics  for  which  I   >   4.   Therefore,  if  I  >  4,    the  integral 
vanishes  by  the  orthogonality  of  the  spherical  harmonics.   By  a  similar 
argument,  all  terms  for  which  H   is  odd  have  zero  matrix  elements  because 
the  density  u*v  has  inversion  symmetry,  and  the  term  U  in  the  potential 
expansion  does  not  have  this  symmetry  for  odd  I,    so  that  the  above  inte- 
grand is  an  odd  function  which  vanishes  upon  integration. 

For  convenience  in  the  subsequent  treatment  of  the  trigonal  distor- 
tion, the  polar  axis  for  determining  the  form  of  the  crystal  field  per- 
turbation is  taken  along  the  [ill]  direction.   This  direction  is  a  three- 
fold symmetry  axis  of  the  octahedron,  and  the  potential  expansion  for  a 

1  Q 

regular  octahedron  in  this  case    is 


V.    =   D'|r4   y0    (8,    t)    +    (^   [y3    (6,    *)    -y"3    (9,    4»H 


(7) 


where  D'  is  the  only  A  coefficient  left  in  the  expansion  (6)  because  of 
the  high  symmetry.   The  elements  in  the  secular  determinant  are  the  in- 
tegrals 

■/>"*  |E  -eV.(r,e,<j>)J  *J'  di,  (8) 

where  the  sum  is  taken  over  all  electrons  outside  closed  shells,  i.e., 
the  potential  function  is  a  single  electron  operator. 
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II>B< Determination  of  Correct  Zero  Order  Wave  Functions 

To  evaluate  the  integrals  (8),  it  is  necessary  to  determine  the 
correct  zero  order  wave  functions  which  are  eigenstates  of  the  free  ion 
Hamiltonian,  not  including  spin-orbit  coupling.   They  are  linear  combi- 
nations of  the  antisymmetrized  single  electron  wave  functions.   In  deter- 
mining these  functions,  use  is  made  of  the  theorem  in  quantum  mechanics' 
which  states  that  if  an  operator  A  commutes  with  the  Hamiltonian  with 
eigenstates  <^,  then  A^  is  also  an  eigenstate  of  the  Hamiltonian.19   m 
particular,  the  lowering  operator 

L_  =  L  -  iL  , 
x     y 

where  L  and  L  are  the  x  and  y  components  of  the  angular  momentum  oper- 
x      y 

ator,  commutes  with  both  the  basic  part  of  the  atomic  Hamiltonian  and 

20 
the  repulsion  perturbation  part.    If  a  wave  function  is  chosen  with 

maximum  magnetic  quantum  number  m  consistent  with  the  orbital  angular 

momentum  quantum  number  SL,    then  all  2£+l  wave  functions  characterized  by 

|m|  <_  I   can  be  generated  by  successive  applications  of  L_.   For  three  d 

electrons,  a  proper  F  state  eigenfunction  having  I   =  m  =  3  is 

+++ 
|£m>  =   1 33>  =  (210),  (9) 

where  the  numbers  in  parentheses  indicate  the  m  values  for  each  of  the 
three  electrons  and  the  +  signs  indicate  that  the  spin  of  each  electron 
is  taken  to  be  +1/2.   This  choice  is  made  for  ease  of  computation  of 
the  crystal  field  energy  levels,  and  since  the  crystal  field  operator 
does  not  affect  spin,  no  loss  of  generality  occurs  if  the  spins  are 
chosen  in  this  way.   These  signs  may  be  omitted  in  later  references  to 
wave  functions  of  this  type,  so  that  electronic  m  values  shown  in  paren- 
theses may  also  be  taken  to  mean  that  the  spin  quantum  number  m  =  +1/2, 
unless  otherwise  noted. 
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The  required  antisymmetric  nature  of  the  wave  functions  is  implicit  in 
the  notation  (210) ,  because  this  expression  is  shorthand  for  the  linear  com- 

K-    *'    21 

bination 

7=  rX2(l)  Xl(2)  XQ(3)  +  X2(3)  Xl(l)  XQ(2) 

+  X2(2)  Xl(3)  x0(l>]   "   ^y[x2(2)  Xl(D  Xq(3) 

+  X2(D  Xx(3)  x0(2)  +  X2(3)  Xx<2)  X0(l)l  . 

in  which  the  subscripts  represent  angular  momentum  m-values,  and  the  numbers 
in  parentheses  are  electron  coordinates.   Note  that  an  odd  number  of  co- 
ordinate permutations  changes  the  sign  of  this  wave  function  and  an  even 
number  leaves  the  sign  unchanged,  as  required  of  an  antisymmetrized  wave 
function. 

m 
The  lowering  operator  L_,  acting  upon  a  wave  function  Xp   gives 

r  11/2 

L_  x£m  =  *   iU+1)  -m  (m-1)      x™      .  (10) 

The  operator  L_  is  defined  as  L  -  +  L,-  +  L  - ,  where  the  subscripts  refer 
to  the  three  electrons  outside  closed  shells.   The  lowering  operator  is 
usually  defined  as  L_/n,  and  the  factor  ft  will  be  omitted  for  conciseness 
in  future  references  to  operators  of  this  type.   The  effect  of  this  operator 
on  the  wave  function  (9 )  is  then 

|32>  =  L_(210)  =  /lO(llO)  +  /l2(200)  +  /l2~(21-l)  .         (11) 
The  first  two  terms  on  the  right  hand  side  of  this  equation  are  forbidden 
by  the  Pauli  exclusion  principle  with  the  result  that  the  m  =  2  wave 
function,  after  normalization,  is 
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|3^>  =  (21-1)  . 
The  other  4F  eigenstates  of  the  repulsion  interaction,  which  are  the 
proper  starting  functions  for  the  crystal  field  calculation,  are  ob- 
tained by  successive  applications  of  the  lowering  operator  L_.   These 
functions  are: 

1 31>    =    J^    [2(21-2)  +  /5  (20-1)]  , 

|30>    =  y|   [(10-1)  +  2(20-2)]  , 

>.       i   P  n  (12) 

|3-J>   =  fi£    [2(2-1-2)  +  /6  (10-2)J  , 

|3-2^>   =  (1-1-2), 

|3-3^>   =  (0-1-2)  . 
The  P  eigenstates  of  the  repulsion  interaction  will  be  required  for  the 
determination  of  the  matrix  elements  of  the  crystal  field  between  4F  and 

4P  multiplets.   These  eigenstates  have  been  determined  by  Theissing  and 

„   -1   20 

Caplan   and  are  given  by 

|ll>  =  ^q  [i/6(21-2)  -  2(20-1)]  , 

|10>  =  ^      [-2(10-1)  +  (20-2)]  ,  (12a) 

|l-£>=  J^   [^6 (2-1-2)  -  2(10-2)]  . 

4 
The  matrix  elements  of  the  crystal  field  perturbation  (7)  between  P 

states  are  all  zero. 
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4 
II. C.   Calculation  of  F  Crystal  Field  Matrix  Elements 

4 
The  7x7  F  crystal  field  perturbation  matrix  is  formed  by  using 

the  wave  functions  (12)  as  basis  functions  for  the  calculation  of  the 

matrix  elements  of  the  perturbation  potential  (7) .   Let  y  be  a  single 

m 

electron  wave  function  with  magnetic  quantum  number  m.   Since  the  per- 
turbation potential  is  a  single  electron  operator,  it  affects  only  one 
of  the  wave  functions  in  parentheses  at  a  time.   For  example,  the  3,3 
matrix  element  of  the  crystal  field  perturbation  is 
<33|v|33>  =  /(210)*V(210)dx 

=  /X2*X2dT  /x1*X1dT  /x0*Vx0dT 

(13) 
+  J*X2*X2dT  /x]_*Vx1dT  fx0*XQ^ 

*   /x2*vX2dx  /x1*X1dx  /x0*X0dx. 

The  integrals  in  which  the  potential  V  does  not  appear  are  all  unity  be- 
cause the  wave  functions  are  orthonormal,  with  the  result 

<J3  Jv|  3^>  =  /x2*Vx2dx  +  /Xl*VXldT  +  /X()*Vx0dT.  (14) 

This  property  of  the  potential  aids  in  reducing  the  number  of  integrals 
which  must  be  evaluated  in  determining  the  off-diagonal  matrix  elements 
as  follows.   Consider  the  integrals 

/xa*xbdT  /xc*XddT  /Xe*VxfdT,  (15) 

and  assume  that  fy   *Vx,-dT  is  non-zero.   The  term  (15)  will  then  vanish 

e   r 

unless  a  =  b  and  c  =  d,  which  implies  that  matrix  element  terms  of  the 
form 

/(ace)V(bdf)dt, 
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which  occur  only  in  off -diagonal  elements,  will  be  zero  unless  at  least 
two  sets  of  functions,  such  as  a,  b  and  c,  d  are  equal.   The  other  ele- 
ments of  the  crystal  field  matrix  are  calculated  as  indicated  in  Equa- 
tion  (13)  . 

Since  V  and  x  are  functions  of  r,  6,  and  <j>,  the  integrands  of  the 

f  (rj  r  where  f  (r) 
is  the  product  of  the  radial  dependences  of  the  wave  functions  x»  which 
are  in  turn  products  of  these  functions  of  r  and  the  y  .(8  #<(>).   Since  the 
exact  form  of  these  functions  of  r  is  not  known,  it  is  conventional  to 

take  the  integrals  /f(r)r  dx  =  rl   as  part  of  the  crystal  field  parameters 

m~F 
to  be  determined.   These  parameters  are  then  of  the  form  A^  r  .   In  mak- 
ing this  assignment,  the  assumption  is  made  that  the  functions  f(r)  are 
constant  over  the  range  of  energies  spanned  by  the  3d  configuration, 
which  is  about  76,500  cm"  .    This  assumption  will  be  valid  as  long  as 
the  3d  configuration  does  not  interact  appreciably  with  higher  energy 

configurations,  an  approximation  which  will  be  accurate  if  the  energy 

3 
differences  between  the  3d  and  higher  configurations  is  sufficiently 

large.   The  magnitudes  of  these  energy  differences  are  not  known,  so  that 

the  accuracy  of  the  above  approximation  must  be  tested  by  agreement  of 

measured  and  calculated  energy  levels. 

A  method  of  general  applicability  for  determining  crystal  field  matrix 
elements  involves  the  use  of  Wigner  coefficients.   In  calculating  matrix 
elements  by  this  method,  use  is  made  of  extensive  tables  of  the  Wigner  3-j 
symbols  available  in  the  literature.   The  matrix  elements  of  the  crystal 

field  perturbation  (7)  were  determined  by  using  Wigner  coefficients  in  the 

22 

following  way.  '   The  matrix  elements  of  the  tensor  operator  T(kq)  between 

state  vectors  of  a  system  which  are  simultaneous  eigenvectors  of  the  angular 

2 
momentum  operators  J   and  J  of  the  system  are  given  by 
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(Y,J'm'|T(kq)|Yjm)  =  (-1)3  _m  /  j 'kj\  (y1  j  '  |  |t  (k)  |  |  Yj)  ,     (16) 

y-m '  qm J 

where  y   is  a  general  quantum  number,  (  j 'kl\     is  a  Wigner  3-j  symbol  and 

V -m  qm/ 

the  final  term  is  the  reduced  matrix  element.   The  spherical  harmonic  oper- 
ator Y£(0,«())  appearing  in  the  crystal  field  expansion  is  a  tensor  operator 
with  k  =  £  and  q  =  m,  as  shown  in  Reference  22.   Since  j  =  I   for  crystal 
field  calculations  for  the  chromium  ion,  the  reduced  matrix  elements  are 
given  by 


U'||y<k)|lA)  =  (-I)1  f"(2^+l)(2fcfl)(2£+l)]^  /Tk£\ 

L     4tt  \o  oq|  ' 


(17) 


where  Y(k)  =  Y°. 


The  definition  of  the  3-j  symbols  provides  another  method  of  re- 
ducing the  number  of  integrals  which  must  be  evaluated  in  determining 

the  crystal  field  matrix  elements.   The  defining  equation  for  the  3-j 

22 

symbols  '  contains  a  term  6(m'-m,q)  so  that  all  of  the  integrals  (15) 

for  which  m  ^  m'-q  are  zero. 

As  examples  of  the  use  of  Equations  (16)  and  (17)  for  calculating 
matrix  elements,  both  a  diagonal  and  an  off-diagonal  element  of  the  crystal 
field  matrix  will  be  calculated.   Diagonal  elements  are  given  by  the  Y 

term  in  the  expansion  (7)  and  off -diagonal  elements  are  given  by  the 

.   3     -3 

<Y4   -  Y4   )  term.   Because  of  the  6(m'-m,q)  term  in  the  definition  of 

the  3-j  symbols,  the  off  diagonal  term  will  couple  only  wave  functions 

which  differ  in  m  values  by  ±3.   Using  Equations  (16)  and  (17),  the  3,3 

matrix  element  is  determined  to  be 
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<33|V4°|33>  =  (2|,Y(4)||2)[(-)-(-)+(^J|-      US, 

The  double-bar  term  (2| |y(4) | | 2)  is  common  to  both  diagonal  and  off- 
diagonal  elements ,  and  is  given  by 

r  -i1/2 

t2||Tc«ii2,-paiaja.]    (-)   . 

Making  use  of  the  table  of  3-j  symbols  in  Reference  22  to  evaluate 
the  3-j  symbol  on  the  right-hand  side  gives 

1/2 


^3157ry 


(2||Y(4)||2)  =f- 


Substitution  of  this  result  into  Equation  (18) ,  and  evaluation  of  the 
3-j  symbols  inside  brackets  finally  gives  for  the  matrix  element 


<»">*> -^PW 


(21) 


18 
The  coefficient  D'  in  the  expansion  (9)  is  given  by  Hutchings   as 

D.  . .  4*  £a  r  (22) 

a 

where  Z  is  the  atomic  number,  e  is  electronic  charge,  a  is  the  distance 

4 

of  the  octahedral  charges  from  the  chromium  ion,  and  r  is  the  average 

value  of  the  fourth  power  of  the  distance  from  the  chromium  ion  at  which 

the  potential  is  measured  and  is  given  by 

~4~         4 

r  =  /f(r)r  dx ,  (23) 

as  mentioned  earlier.   Muliplying  Equation  (21)  by  D'  gives  the  matrix 

element 

-  {?  <£)  (ill)    <2>  ■  -»* 
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The  first  and  second  terms  in  parentheses  in  Equation  (24)  are  conven- 
tionally designated  D  and  q  respectively,  and  are  the  basis  for  specify- 
ing energy  levels  in  all  iron  group  ions.   The  constant  Dq  is  usually 
determined  from  optical  absorption  measurements  as  will  be  discussed  later. 

As  an  example  of  the  calculation  of  an  off -diagonal  matrix  element, 
the  term  3,0  will  be  evaluated.   Making  use  of  the  wave  functions  of 
Equation  12,  this  matrix  element  is  given  by 

<J3|  (Y43  -  Y4~3)  |30>  =  ■=  /(210)  (Y43  -  Y4"3)  (10-1)  dx 

+  £  /(210) (Y  3  -  Y  ~3) (20-2)dx.  (25> 

/5  4      4 

Note  that  in  the  first  integral,  m1  =  2  and  m  =  -1 ,  while  q  may  be  either 

-3 
±3.   The  requirement  stated  earlier  that  m  =  m'-q  will  then  cause  the  Y4 

term  in  the  first  integral  (27)  ,  to  vanish.   In  the  same  way,  it  is  easily 

-3 
seen  that  the  Y    term  in  the  second  integral  also  vanishes.   In  computing 

the  0,3  matrix  element,  it  is  found  that  the  Y   term  vanishes  in  each  inte- 
gral, but  this  element  is  equal  to  the  3,0  element  in  accordance  with  the 
requirement  that  quantum  mechanical  operators  are  Hermitian. 

In  computing  the  off-diagonal  elements,  Equation  (16)  is  used  in  the 
form 

(yj 'm' |T(kq) | yjm) 

=  (-l)J-m  <3'J»,-»lj'jP"i>  (j'||T(k)||j),  (26) 

(2k+l)  ' 


in  which 


(j'jm'-mlj'jkq)  =  (-l)j'~3  q(2k+l)  [  J',kj  ) 

V-m  qm J 
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is  called  the  vector-coupling  coefficient.   Note  that  the  reduced  matrix 
element  in  Equation  (26)  is  the  same  as  that  of  Equation  (16).   The  3,0 
element  is  then 

<"|V43|30>=l(^)(3k)'![-/|(2221|2243, 

-   1   (2212(2243) I  .  (28) 

/5  -1 

Using  the  tables  of  vector-coupling  coefficients  found  in  Condon  and 
Shortley  "  and  multiplying  by  L~-J  D'  as  required  by  Equation  (7) 

finally  gives 

<J3  |v|  30>  =  2  /lO~Dq.  (29) 

The  other  matrix  elements  are  calculated  in  a  similar  way,  so  that  the 

4 
7x7   F  crystal  field  perturbation  matrix  is  found  to  be 

m£       0       3       -3       1       -2        -12 

0  -4  2/l0        -2/To 
3             2/l0        -2  0 

-3  -2/LO  0  -2  (3°) 

1  -2/3     -4/3/5 
_2                          -4/3/5   14/3 

-1  -2/3      4/3/5 

2  4/3/5     14/3 
where  all  of  the  elements  are  expressed  in  terms  of  the  constant  Dq. 

Solution  of  the  3x3  submatrix  gives  E  =  -12Dq,  -2Dq  and  6Dq.   The 
other  two  submatrices  each  give  E  =  -2Dq  and  6Dq  so  that  the  4F  energy 
levels  for  the  chromium  ion  in  a  purely  octahedral  field  are  found  to  be 
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E   =  -12Dq  singlet, 

E   =  -2Dq  threefold  degenerate,  (31) 

E   =   6Dq  threefold  degenerate. 

The  state  functions  for  each  of  the  levels  (31)  are  found  by  sub- 
stituting the  eigenvalues  into  the  original  matrix  equation,  solving  for 
the  coefficients  of  each  of  the  component  states,  and  normalizing  these 

coefficients.   For  example,  the  coefficients  C  of  the  r_  state  are 

v  m         2 


(32) 


giving  the  three  equations , 


-4CQ  +  2/To  C3  -2/I0  C_3  =  -12  C  , 

2^L0  CQ  -2  C3  =  -12  C3, 
-2/l0   CQ  -  2C_3  =  -12  C_3. 


Solving  for  the  coefficients  and  normalizing  gives,  for  E  =  -12Dq, 

jS,;?s  _^\  .  .5. 

3 


>  4|30>4|33>+4|3-3>. 


(33) 


where  the  Roman  numeral  is  a  designation  for  keeping  track  of  the  various 
functions.   In  a  similar  way,  the  other  states  are  determined  to  be: 
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E    =    -2Dq 
">    -   ||33>   ♦   /i|3-3>    , 

"I>-/||3>     +/||3-2>, 
>    "   /f|3-l>    -    4l30    - 


IV 


E    =   6Dq 


(33a) 


V>      -||30>    +i     /I   (|33>   -    |3-3>    ), 
VI>     =/i|31>    -      /|j3-2>, 
VII>-/I|3-1>    +   /||32>. 


In  Equations  (33a),  the  first  three  functions  are  the   F   states,  and 

4 
the  last  three  are  the  F   states.   In  the  expressions  for  the  6  Dq  wave 

functions,  the  primes  are  used  to  distinguish  these  wave  functions  from 

4 
those  that  will  result  when  the   P  interaction  is  considered.   These  wave 

functions  are  the  U  representation  for  the  F  multiplet.   Since  the 

crystal  field  perturbation  has  no  effect  on  the   P  multiplet  energy  level, 

Equations  (12a)  give  the  Lr  representation  for  that  multiplet.   The  wave 

functions  (12a),  (33)  and  (33a)  will  be  used  as  a  basis  for  evaluating  the 

4       4 
interaction  of  the  F  and  P  states  through  the  crystal  field  perturbation. 

4        4 
II. D.   Interaction  of  the  FT   and  P  States 

3+  4  -1 

In  the  free  Cr   ion,  the   P  level  lies  about  14,000  cm   above  the 

4  -1 

F  ground  level,  but  this  interval  is  reduced  to  about  10,000  cm   when  the 

24 

ion  is  placed  in  a  crystalline  field.   Owen   has  shown  that  this  reduction 
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in  energy  is  caused  by  covalent  bonding  between  the  Cr  +  ion  and  the 
surrounding  oxygen  atoms,  reducing  the  coulomb  interaction  between  the 

d  electrons  and  resulting  in  a  reduction  in  the  apparent  energy  level 

4 
of  the   P  state.   Covalent  bonding  has  also  been  considered  by  Sugano 

25 
and  Peter   in  their  treatment  of  the  energy  spectrum  of  ruby.   The 

4 
crystal  field  perturbation  (7)  has  no  effect  on  the   P  energy  level  be- 
cause the  matrix  elements  of  this  perturbation  in  the   P  state  are  zero. 
However,  after  the  application  of  the  perturbation  (7)  to  the  ground 
level  and  the  determination  of  the  eigenfunctions  (33a)  ,  it  is  found 
that  the  crystal  field  perturbation  couples  the  FT   states  and  the   P 
level  with  resultant  changes  in  energy  levels  and  mixing  of  states. 

4 
The   P  wave  functions  for  the  free  chromium  ion  are  given  by  Equations 

(12a)  .   Each  of  these  functions  is  coupled  to  one  of  the  three  FT.  wave 

functions  by  the  perturbation  (7) .   The  matrix  elements  of  V  between  each 

4  4 

of  the  three   P  states  and  the  corresponding  FT   state  are  determined  in 

the  same  way  as  were  the  elements  for  the  matrix  (30) ,  and  are  found  to 

be 

<Cy'    |v|l0>  =  4  Dq, 

<VI'   |vjl£>  =  <VII'  |v|l-l>  =  -4  Dq,  (34) 

in  the  notation  of  Equations  (12a)  and  (33a) .   These  results  differ 
slightly  from  the  results  of  Finkelstein  and  Van  Vleck  ,  who  obtained 

4  Dq  for  each  of  the  matrix  elements  (34) .   Note  that  this  sign  difference 

4      4 
will  not  affect  the  energy  interval  between  the  admixed  F  and  P  levels, 

but  it  will  affect  the  state  compositions  of  these  levels.   The  matrix 

4      4 
which  accounts  for  the  coupling  between  P  and  FT   levels  is  then,  in 

multiples  of  Dq, 
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H 

H 

4-4 

H 

4-2 

-12 

4-s 

-2 

H 

6 

±4 

X 

±4 

E(4P) 

(35) 


where  E(4P)  is  the  energy  of  the  4P  state  of  the  ion  measured  from  the 
ground  state  (i.e.,  from  4F  =  0) .   Solution  of  this  matrix  yields  a  sing- 

A  A  A 

let  Fr2  at  -12Dq  and  a  triplet  F5  at  -2Dq  as  before,  but  the   FF4  and 

4 
Pr4  levels  are  shifted  by  the  perturbation.   The  positions  E4  and  E5  of 

these  new  levels  are  found  by  solving  the  2x2  submatrix  in  Equation 
(35)  and  are 

1*5 


[(Effil^ica)  '  +  (4  ^ 


.  6  Dq  +  E(P)     fE(P)  ~  6  Dg 
E4  ~      2 V 2 )       +  (*   Uq;  I!  '  (36) 


and 


_  6  Dq  +  E(P)    HE_(P)_-6Dq\  2  J   .  .  ^  N2 

E5 2 \\ 2 J      (4  Dq)    *  (37) 

Figure  2  shows  the  effects  on  the   F  and  P  levels  of  the  Cr   ion  which 
result  from  the  application  of  the  perturbation  (7) .   The  energies  of  the 
three  excited  levels  shown  in  this  figure  may  be  determined  from  optical 
absorption  measurements  and  the  constants  Dq  and  E(P)  may  be  calculated 
based  on  these  measurements . 
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4      r 

E  =  E(P)      P        4  /     Energies  Given  by 


f-       E  =  6Dq  X  Equations  (36)  and  (37) 


4 
E  =  0    F 


E  =  -2Dq 


T  E  =  -12Dq 


4      4 

Figure  2.    F  and  P  Energy  Levels 

of  Cr3+  in  YAG 

It  is  convenient  at  this  point  to  look  ahead  to  the  experimental 
results  for  a  solution  to  Equation  (36)  for  E(P),  because  the  value  of 
this  parameter  will  be  such  that  an  exact  solution  of  the  2x2  sub- 
matrix  of  (35)  will  be  obtained.   In  Section  V,  it  will  be  shown  that 

4  -1  -1 

the  FT   and  E  levels  lie  at  16,700  cm   and  23,300  cm   respectively. 

4 
The  location  of  the  FT   level,  which  is  lODq  above  the  ground  level, 

gives  Dq  =  1670  cm   .   Equations  (31)  show  that  the  free  ion  level  lies 

12Dq  =  20,000  cm   above  the  ground  level.   The  E   level  then  lies 

-1  4 

23,300  -  20,000  =  3300  cm   above  the  free  ion  F  level  and  this  sepa- 
ration is  equal  to  2  Dq  within  experimental  error.   To  summarize  the 

above  discussion,  the  energies  of  the  pertinent  levels  relative  to  the 

4 
free  ion  F  level  are 
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Fr2:   -12Dq  =  -20,000  cm  1 

4  -1 

Fr5:   -  2Dq  =  -3340  cm  (38) 

E4   :     2Dq  =  3340  cm"1. 

Equation  (36)  may  be  solved  for  E(P)  to  give 

(E   -  8Dq)  (E   +  2Dq) 

E(P)  =  — % .  (39) 

E4  -  6Dq  (Jy) 

Substituting  E^  =  2Dq  into  this  equation  gives  E(P)  =  6Dq  =  10,000  cm"1, 
in  good  agreement  with  similar  calculations  performed  on  chromium  in  three 
different  complexes  as  cited  by  Low,   who  gives  F(P)  =  10,000  cm"1  for 
Mg0:Cr   and  CrfNH^  +,  and  10,200  cm"1  for  Cr(H20)  3+.   when  substituted 
into  Equation  (37),  the  above  values  of  Dq  and  E(P)  give  E  =  lODq  = 

16,700  cm   .   The  calculated  energy  levels  for  Cr  +  in  a  purely  octahedral 

4   4 
field,  including  F/  P  interaction,  are  then 

E    =  -12Dq  singlet  (0) , 
2 

E„   =  -2Dq  threefold  degenerate  (16,700), 

5 

E    =  2Dq  threefold  degenerate  (23,400),  (40) 

41 

E    =  lODq  threefold  degenerate  (36,700), 

42 

where  E    and  E    denote  the  energies  of  the  lower  and  upper  r   levels 
41       42  4 

respectively.   The  numbers  in  parentheses  are  the  energies  in  cm"  relative 
to  the  ground  level,  and  will  be  compared  to  experimental  values  in 
Section  V. 
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The  state  functions  for  each  of  the  energy  levels  (40)  are  found  as 
before  by  substituting  the  eigenvalues  into  the  original  matrix  equation, 
solving  for  the  coefficients  of  each  of  the  component  states,  and  normal- 
izing these  coefficients  to  unity.   The  T     and  r   states  are  not  affected 

by  the  perturbation  V  because  matrix  elements  of  V  between  these  states 

4 
and  the   p  states  are  zero.   Therefore,  it  is  only  necessary  to  use  the 

4   4 
2x2  submatrix  of  (35)  to  determine  the  compositions  of  the  admixed  F/  P 

levels.   In  making  this  calculation,  the  sign  of  the  off-diagonal  element 


is  chosen  to  be  consistent  with  Equation  (34] 
wave  function  |v^  ,  the  matrix  equation  is 


For  the  2Dq  level,  and  the 


~V2 


p2 


=  2 


Vw 


(41) 


where  C   and  C   are  the  coefficients  of  the  component  states  for  2Dq. 


This  set  of  equations  gives  C. 


-  -n        =     1   ,  so  that  the  new  wave  function 
V2     p2    — 


is 


V 


>  =  J_|V>  -Jh io> 

/2  Si 


=   ^|30>+^|33>-f 


3-3>    -   -^|10>     . 
•2 


(42) 


In  a  similar  way,  the  remaining  wave  functions  for  E  =  2Dq  (T   )  are  found 
to  be 


|VI>      --^|31>-i     /||3-2>+^|U>     , 

2/1  /2 


1>.^|3-1>+i/||32>      +-i|ll-l>    , 

2/1  *        J  /2 


(42a) 


29 


and  the  states  for  E  =  lODq  (r   )  are 


.  ^U«\  .  /5,_\    •§-, 


IX 


1>=^3»>^|33>-f|3-3>t-i|10>, 

V2 

These  states,  along  with  the  T2  and  Tg  states  given  by  (33)  and  (33a)  will 
be  used  to  determine  the  effects  of  the  trigonal  distortion  on  each  of  the 
threefold  degenerate  energy  levels. 

II. E.   The  Effect  of  the  Trigonal  Distortion  in  First  Order 

The  trigonal  distortion  of  the  crystal  field  in  YAG:Cr,  resulting  from 
a  slight  displacement  of  the  chromium  ion  along  the  crystalline  fill]  direc- 
tion, has  a  considerable  effect  on  the  location  of  the  crystal  field  energy 
levels.   This  distortion  manifests  itself  as  a  splitting  of  a  few  hundred 
wave  numbers  in  the  energies  of  the  r  and  r   states.   Since  these  levels 
lie  16,000  to  20,000  wave  numbers  above  the  T     ground  level,  it  is  convenient 
to  consider  this  trigonal  distortion  as  a  perturbation  of  the  basic  states. 

The  crystal  field  expansion  in  terms  of  spherical  harmonics  which  in- 
cludes the  effects  of  the  trigonal  distortion,  may  be  written11 

,0   2   0     0   4   0     3   4    3     -3 
V  =  A2  r  y2  +  A4  r  Y4  +  A4  r  {y^  -  y^   ) .  (44) 

For  a  crystal  field  of  octahedral  symmetry,  the  potential  (44)  will  approach 
the  basic  potential  (7)  as  the  coefficient  of  the  y   term  approaches  zero. 
It  will  be  shown  in  Section  V  that  the  calculation  of  crystal  field  energy 
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levels  using  (7)  gives  results  that  agree  within  about  3.9  percent  with 
the  measured  centers  of  gravity  of  the  optical  absorptions.   These  results 
imply  that  the  effects  of  the  trigonal  distortion  on  the  locations  of  the 
centers  of  gravity  of  the  various  energy  levels  is  negligible,  and  that  the 

contribution  of  the  y   term  is  small.   If  this  contribution  is  small,  the 

0       3-3 
coefficients  of  the  y  and  (y   -  y   )  terms  are  almost  unaffected  by  the 

trigonal  distortion.   The  crystal  field  potential  expansion,  including  the 

trigonal  distortion,  is  then  given  to  good  approximation  by 


0   2   0     ,4 
V  =  A„  r   y„  +  D'  r 
2    ^2 


0    ,lo//2  ,  3    -3, 

y4  +  (->         (^4  "  y4  > 


(45) 


0   2   0 
The  term  A  r   y   is  then  treated  as  a  perturbation  of  the  basic  states  (33) , 

(42)  ,  (42a)  ,  and  (43)  ,  which  are  eigenstates  of  the  basic  potential  (7) . 

The  states  T    ,    T      ,  T.      and  r  are  separated  by  energies  on  the  order 
of  10,000  cm   ,  and  it  is  therefore  feasible  to  consider  the  trigonal  dis- 
tortion as  a  perturbation  of  each  of  these  states  individually.   To  account 
for  interaction  between  states,  it  is  necessary  to  carry  the  perturbation 
calculation  to  second  order.   The  r  wave  functions  interact  with  the  T 
ground  level  wave  functions  through  second  order  spin-orbit  coupling, 
splitting  the  F_  level  into  two  sublevels  separated  by  less  than  one  cm 
This  splitting  can  be  measured  by  the  method  of  PMR,  and  if  the  spin-orbit 
coupling  coefficient  is  known,  the  splitting  of  the  T      level  due  to  the 

trigonal  distortion  can  be  calculated.   The  value  of  this  splitting  is  then 

0   2 
used  to  evaluate  the  coefficient  of  the  trigonal  distortion  A  r  .   Using 

this  result,  the  effects  of  the  trigonal  distortion  on  the  V     and  T      levels 

can  be  determined . 
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The  matrix  elements   of  the   trigonal  distortion 

0      2      0 
Vtrig  =  A2   r     ^2  (46) 

will  now  be  evaluated  using  the  F  wave  functions  as  the  basis.   This  cal- 
culation is  made  in  the  same  way  as  the  matrix  elements  of  the  basic  poten- 
tial were  calculated,  and  the  desired  matrix  elements  are  determined  to  be 

O^trigl1^  =K33lVtrigl33> 

+  ^<3-3|v   .  |3-3>  -i-  (^\\°7  , 
2      >.   i  trig1  *<        14  \t\J        2 

<$n|v.    .   |iii>=<jv|v        |iy>=-<3i|v        |31> 

>»        '    trig1      ^      ^s     l    trig1     ^        6  V    '    trig1    ^ 


(47) 


+  K3-2ivigi3-2> +  ¥  <3ii  vigi3-2> 

,5//2  ,1,02 

=  -(7}    WA2r- 

Note  that  all  off-diagonal  elements  are  zero  because  y   is  a  diagonal 
operator.   Equation  (47)  therefore  gives  the  energies  of  the  singly  and 
doubly  degenerate  sublevels  of  the  T      state,  correct  to  first  order  of  the 
trigonal  distortion. 

The  trigonal  distortion  also  splits  each  of  the  two  V      levels  into 

doubly  degenerate  and  nondegenerate  levels  in  the  same  way  as  the  r   level 

5 

is  split.   The  first  order  energies  of  the  component  states  of  the  r   level 
with  E  =  2Dq  (V      )  are  determined  to  be 


32 


<V|*»   r2  y°|v>   .  4   r2    [f  <30|y°|30>    +   5-  <33|y» 

*  h  <«!*>£>   ♦  \  <10|y°H0>    -  f  <30|y»] 

31     f  jA  "*      0  ~~2 
"        140   ^ttJ        A2   r      ' 

(48a) 

<;      |0      20i\  s*        I    0      2      0 1       \. 

VI I A     r     y    |vi>    =  <Vll|A     r     y    |vil> 


0      2 
-A2r 


Qj<hl»»+!2<»hr» 


Kl-1|yo|l.1>+^<3l|yo|l>]=.i_^\o 


In  a  similar  way,  the  energies  of  the  T      level  with  E  =  lODq  (V      J  are  given 
by 


<,  0   2   0 1    V      1  ^5V  0   2 
VIIl|A2r  y2|vill>  =-^[-j     A2r   , 

<'  i  0   2  0i  \    /|0   2   0|\    1  YsY5   0 
IX|A2  r  y2|lX>  =  <X|A2  r  y2£>=  —  ^J  A2 


(48b) 


*-2 
r 


The  energy  levels  of  all  of  the  T     and  T     states,  correct  to  first  order 

0   2   0 
of  the  perturbation  A  r  y  ,  are  then 


En  r2 


14  \^J   A2 

28  ^j    2 

(f) 

31    (s\  0   2 

=  2D*  "  280   V~J   A2  r   ' 


=  EIV  =  "  I   A°  ^ 


r  5  \  ^    0   2 

2Dq  +  nb"   ?      A2 r 


E      =  E 
VI       VII 


(49) 


ii 


eviii 


!80  \k)        A2 


Eix    =  Ex  =  lODq  +  2=5-1   j 


II. F.   The  Effect  of  the  Trigonal  Distortion  in  Second  Order 
The  T5  and  the  lower  I*  states  are  spaced  4Dq  =  6680  cm"   apart,  so 
that  it  is  reasonable  to  expect  that  second  order  trigonal  distortion 

effects  might  influence  the  positions  of  the  energy  levels  and  affect 

0  ~~2 
determination  of  the  coefficient  A  r  .   The  second  order  correction 

E  '   to  the  energy  E  of  a  state  |r^  caused  by  a  perturbation  V  is 

p(2)       y    I  <n|v[k>  I2 

,  ,  E,  -  E      '  (^U) 

Mn      k    n 

where  the  index  k  refers  to  other  states  coupled  to  |  n^>  by  the  pertur- 
bation.  To  apply  Equation  (50)  to  the  V      states,  it  is  necessary  to  de- 
termine matrix  elements  of  y  between  the  r   states  and  all  of  the  other 
states.   A  careful  examination  of  the  eigenfunctions  (33),  (42)  and  (43) 
shows  that  many  of  these  matrix  elements  are  zero.   For  example,  the  matrix 
element  <^l|y  l^is  given  by 

-i<5^>>+i<M^iH>- 

But   <Q3|y2|33>    =  \3-3|y2|3-3>    so   that  the  abOVe  matrix  element   is 
zero.      Similarly,    <Il|y2|y>    and    <Il|y2|viII>    are  also  zero.      The  only 
nonvanishing  matrix  elements   are  determined  to  be 
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and 


where 


<m|y°|vi>  -  <iv|y°|vn>  .  i  /|  <3i|y»|3> 
4/f<Hy°|3-2>-i/!<3i|y°|n> 

20/2     \t\J     ' 
<m|y°|ix>  -  <iv|y°|x>  .  1  /f  <3l|y°|3i> 

i/F<H^lH>+t/F<?^°l"> 
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(51) 


(52) 
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These  results  show  that  the  level  of  the  state  1 1^  is  unchanged  by  the 

second  order  correction.   The  energies  of  |lll^  and  |  IV^  are  equal, 
and  are  given  by  the  following  relation,  correct  to  second  order: 


E  =   E 

III  IV 


"@ 


kl) 


2Dq 
2 


28      \ttJ        A2 


160 


31       (5\  h     0 

_4Dq   +    280      [*)        A2 
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\t\J      2-28-1' 
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^02 
A2   r 
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0   2 
The  term  involving  A  r   in  the  first  bracketed  expression  in  Equation 

(53)  will  not  be  small  enough,  compared  to  4Dq,  to  be  negligible,  but 
the  corresponding  term  in  the  second  expression  will  be  negligible  com- 
pared to  12Dq.   The  basis  for  this  assertion  will  be  evident  when  the 

0   2 
final  values  of  Dq  and  A  r  are  calculated. 


The  splitting  A  of  the  r   level,  correct  to  the  second  order  of  the 

0   2   0 
2  r   Y2 

h     n   -=■ 


perturbation  A  r   y  ,  is  then 


A=E    -E    =^[i]   Aur2 
II     III    28 


(f)  *° 


+  ^  (A°r2)   ( 1 _ 

lieo[4Dq+||,(f)       A°r2] 

+ S I 

2- 28-140 (12Dq)  (' 

Substituting  for  the  constants  and  clearing  fractions  gives  the  follow- 

0  ~~ 2 
ing  cubic  equation  for  A  r   in  terms  of  Dq  and  A: 


2.2967  x  10  3 


fo~V  (  0  — \ 

I  A  r   I   +  0.41169Dq  Ia  r   1 


(54) 


+  12Dq  (0.54068Dq  -  0.13967A)  A°  r2  -  48  (Dq) 2  A  =  0. 
This  equation  will  be  solved  in  Section  V  when  values  of  Dq  and  A  are 


determined. 


The  second  order  trigonal  distortion  corrections  to  the  T     and  T 

2  4 

levels  are  determined  in  the  same  way  as  were  the  analogous  corrections 
for  the  T5  level.  For  the  T     level,  the  corrected  energy  is 
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ET    =   -12Dq  - 
I  H        2-7-35."  lirl        14Dq 


7-35  ^Trj 


2  (55) 


2-7-35  VttJ    22Dq 
In  writing  Equation  (55) ,  the  splittings  of  the  levels  to  which  the  r 
level  is  coupled  by  the  perturbation  are  considered  small  compared  to 
the  level  separation.   The  energies  of  the  r  levels  correct  to  second 

order  will  not  be  listed,  but  these  energies  will  be  determined  in  Sec- 

0  ~~ 2 
tion  V  after  the  A  r  coefficient  has  been  calculated. 

Note  that  the  effect  of  the  second  order  y  correction  on  the  V 

ground  state  is  to  slightly  depress  the  ground  state  energy.   This 

change  will  have  a  small  effect  on  the  determination  of  the  constant  Dq 

as  will  be  shown  in  Section  V. 

II. G.   Exact  Solution  of  the  Crystal  Field  Matrix 

Including  4F/^P  Interaction 

and  Trigonal  Distortion 

Although  the  calculations  of  Section  II. F.  give  the  crystal  field 

energy  levels  to  good  accuracy,  it  is  still  meaningful  to  solve  the 

10  x  10  crystal  field  matrix  which  includes  all   F  and  P  states  and  the 

effects  of  the  trigonal  distortion.   This  calculation  was  carried  out 

earlier  in  Section  11. D.  without  including  this  distortion. 

In  making  this  calculation,  the  trigonal  distortion  is  included 
from  the  beginning,  with  the  result  that  all  of  the  resulting  wave  func- 
tions will  be  diagonal  in  this  perturbation.   The  calculation  of  Section 
II. F.  is  therefore  considered  more  accurate  than  the  solutions  of  the 
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10  x  10  matrix  including  the  distortion  because  second  order  effects 
were  shown  to  be  significant  in  the  last  section.   The  most  important 
reason  for  solving  the  10  x  10  secular  equation  which  includes  the  dis- 
tortion is  to  show  how  the  eigenstates  are  mixed  by  this  perturbation. 
In  particular,  it  will  be  shown  that  one  of  the   P  states  is  slightly 
mixed  with  the  ground  state  | ^>  given  by  Equation  (33). 

Again,  several  new  matrix  elements  must  be  evaluated  before  this 
calculation  can  be  made,  namely  the  elements  ^3m  |y Jim /^>.   Thes 
elements  are  evaluated  as  before  and  are  given  by 


ise 


<&!*»  =-#    (f)  \ 


(56) 


The   10  x   10   crystal   field  matrix  breaks  down  into   a  4  x  4   and   two   3x3 
matrices.      The   4x4   matrix   is 

^m.  3,0  3,3 


l 


15V   0   2 


3,0      -4Dq-^^)  AV  2/L0Dq 


35\tt/      2 


3,3  2.fi0Dq  -2Dq4U|)A°r2 


3,-3  -2/lODq 

?q-3-5\7j  "2*  '  3 


i.o  ^-6Ji)\y    +i^b-Dq 


3,-3 

1,0 

-2v/10Dq 

+  3^-35(7)  A2r 

0 

0    2 
A2r 

+  |Vl0Dq 

(57) 

-  rVlobq 

-  |/L0Dq 

-vffi'F  ■ 

The   two   3x3  matrices   are, 


£,m£  3,1 


-2-  ffi  »° 

70   ^ttJ      A2 


,    ,  2  3      /5\        0      2 

3,1  -^Dq  -—  [  a,    r 


3-2  -|/5Dq 


1.1  -  |^Dq  -M  (1)     A?  r2 


(^ 


3-2 

1,1 

-|/5Dq 

-  !**  #(f)'  4  ? 

14,, 

+  2/  ^Oq           (58a) 

+  2/±§Dq 

1      f5V    ,0  ~2 

6Dq  "lb"  7     A2  r    ' 

and 


£,m£  3-1  3   2  1-1 

3-1  -I043  "to  \7j  A2  r  7^Dq  '  ?**■  -S%J    A2  r 

3    2  4/5Dq  i^q  -   2/^Dq  (58b) 

■5  o  j 


1-1 


2  rzr,     2/6  f  5Y  ,0    2      „/io  ,       1    M35    0 

T^  -35-V7J   A2  r    "  2/-l°q       6Dq  -lb"  tJ    A2 


2 

r 


The  results  Dq  =  1670  cm   and  A  r   =  13,300  cm   ,  which  will  be  de- 
rived in  Section  IV,  are  now  used  to  solve  the  secular  Equations  (58) . 
Upon  making  these  substitutions,  the  matrices  (58)  become 


l'\ 

3,0 

3,3 

3,-3 

1,0 

3,0 

-7,639 

10,562 

-10,562 

1,576 

3,3 

10,562 

-2,141 

0.0 

3,521 

3,-3 

-10,562 

0.0 

-2,141 

-3,521 

1,0 

1,576 

3,521 

-3,521 

13,380 

l, 

'mz 

3,1 

3, 

,-2 

1,1 

3, 

.1 

-1, 

,832 

-4, 

,979 

-5,076 

3, 

-2 

-4, 

,979 

7, 

,793 

6,098 

1, 

.1 

-5, 

,076 

6, 

,098 

8,341 

(59a) 


(59b) 
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*'ml 

3,-1 

3,2 

1,-1 

3,-1 

-1,832 

4,979 

-5,076 

3,2 

4,979 

7,793 

-6,098 

1,-1 

-5,076 

-6,098 

8,341 

(59c) 


where  all  of  the  above  matrix  elements  are  expressed  in  cm  .   In  gen- 
eral, only  three  significant  figure  accuracy  is  allowed,  and  the  solu- 
tions to  the  above  secular  determinants  must  be  rounded  off  to  this  ac- 
curacy.  The  accuracy  of  the  measurements  leading  up  to  the  determination 
of  the  elements  of  the  matrices  (59)  will  be  discussed  in  Section  V. 

The  secular  determinants  (59)  were  solved  on  the  Martin  Marietta 

IBM  370  computer  using  a  matrix  diagonalization  program  available  in  the 

program  library.   This  program  calculates  eigenvalues  and  normalized 

eigenvectors  for  complex  matrices  of  dimensionality  up  to  20  x  20.   The 

results  of  this  calculation  are  given  in  Table  I. 

TABLE  I 

Eigenvalues  and  Eigenvectors  of  the  Crystal  Field  Matrix 
Including  the  Trigonal  Distortion 

4   -1 
E  =  -2.02  x  10  cm 

|l'>>        =   0.76l|30>    -    0.457|33>    +    0.457|3-3>    +   0.060|l0> 

E   =   -4.54   x   103    cm-1 
|lll£>    =   0.934|3^>    +   0.255|3-^>    +   0.248|l^ 
J>      =   0.934|3-£>     -   0.255|32>    +   0.248|l-l> 

E   =   -2.14   x   103    cm"1 


IV 


11'^     =   0.707|33>    +   0.707|3-3> 


E   =   6.81   x   10      cm" 


|v'j>       =   0.552|3<£>    +   0.42l|33]>    -    0.42l|3-3>    +    0.584|l£> 

(continued) 
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TABLE  I  (continued) 


E  =  1.96  x  103  cm  1 


VI 


j^>     =  0.011 1  31    -  0.718  |  3-2    +  0.696  1 11 
VII'^>    =0.011 1  3-1    +0.718|32    +0.696|l-l 


|IX'> 

|x-> 


4    -1 
E  =  1.69  x  10   cm 


=  0.355  31    -  0.648  3-2    -  0.674  11 


0.355|3-1    +  0.648|32    -  0.674|l-l 

4    -1 
E  =  1.70  x  10   cm 


JVIII'^>   =  0:342|3<£>  +  0.338|33>  -  0.338|3-3>  -  0.809|l0> 


Comparison  of  the  coefficients  of  the  states  given  in  Table  I  with  those 
derived  without  including  the  trigonal  distortion  shows  that  the  trigonal 

distortion  has  little  effect  on  the  state  compositions.   As  noted  earlier, 

4         i  \ 

there  is  a  very  slight  mixing  of  the   P  state  | 10^>  with  the  ground  state, 

which  will  be  shown  to  have  little  effect  on  the  PMR  results. 


SECTION  III 
THEORY  OF  PARAMAGNETIC  RESONANCE  IN  YAG  DOPED  WITH  CHROMIUM 

III. A.   Perturbation  of  the  Ground  State 
by  the  Magnetic  Field  and  Spin-Orbit  Coupling 

Paramagnetic  resonance  (PMR)  transitions  are  observed  between  the 

m  =  +3/2  and  m  =  +1/2  doublet  sublevels  of  the  I"   ground  state  of  Cr 

s  s  -^ 

in  YAG.   Second  order  spin-orbit  coupling  interaction  between  the  ground 

(T  )  and  the  first  excited  triplet  (T  )  level,  together  with  the  trig- 

onal  distortion,  combine  to  split  the  ground  level  into  two  Kramers 

doublets,  m  =  +3/2  and  m  =  +1/2.   PMR  transitions  are  observed  between 
s   —         s   — 

these  various  doublet  levels,  providing  a  means  of  determining  the  split- 
ting of  the  ground  state  and  the  spin-orbit  coupling  parameter.   These 
multiple  transitions  are  called  the  PMR  fine  structure.   The  results  of 
PMR  measurements  also  give  the  splitting  A  of  the  V     level,  thus  giving 
the  other  condition  for  calculating  the  energy  levels  E  and  E  . 

An  effective  perturbing  Hamiltonian  H   can  be  written  in  terms  of 

the  magnetic  field  H,  spectroscopic  splitting  factor  g  ,  and  spin-orbit 

coupling  parameter  A  as  follows: 

H   =  3H-  (L  +  g  "!")  +  AL-"S",  (60) 

s  s 

where  the  operators  in  this  case  are  state  operators  instead  of  the  single 

electron  operators  that  were  used  in  the  crystal  field  calculations  of 

Section  II.   It  is  convenient  to  evaluate  the  matrix  elements  of  this 

Hamiltonian  using  the  representation  for  L,  with  the  mutually  commuting 

2       2 
operators  L  ,  L  ,  S  ,  and  S  ,  because  both  the  crystal  field  energy  contri- 

2 
bution  and  the  operators  S   and  S   are  diagonal  in  this  representation. 

The  basis  functions  for  this  representation  are  the  eigenfunctions  of  the 

crystal  field  perturbation  with  the  spin  quantum  number  taken  into  account. 

In  particular,  since  PMR  is  observed  only  in  the  ground  state  at  room 
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temperature  without  optical  excitation,  only  the  ground  state  wave 
functions  need  be  considered  in  the  first  order  perturbation  calculation. 
The  four  ground  state  basis  functions,  including  the  spin  quantum  numbers 


S  and  m  are  then 
s 


"C>=  |r2S»>=  4|30|Bs>-4(|33|„s>-  |3-3§n>). 


(61) 


3      1 
In  this  equation,  m  takes  on  the  values  +-  and  +-,    giving  the  four 

S  Z  A 

basis  functions. 

The  matrix  elements  of  the  first  order  perturbation  are 

<T  Sm  |h  |r  Sm  "S  =  <"t  Sm  I  3H  -*L  I  T  Sm  ^> 
•  2   s '  s '  2  s-^        X,  2   s '     '2  s^ 

+  <JnSm  |g  3H-"s|r  Sm  *S  +  A  <*T  Sm  it-sir  Sm  *S  . 
^2   s1  s     '2   s —     V.  2   s1    '2  s^ 

The  components  of  H  may  be  written 

H  =  HcosG, 

z 

+  i  At 
H  =  HsinO  e    , 


(62) 


(63) 


where  6  and  $  are  polar  angles  referred  to  the  symmetry  axis  of  the  crystal, 
which  is  the  |_lll]  direction  for  YAG,  and  is  the  direction  around  which  the 
crystal  field  potential  was  expanded  in  Section  II. A.   Using  the  above  ex- 
pression for  H+,  the  operator  H-L  becomes 

t'"Z  =   H  L   +  i(HL   +  H_L  )  .  (64) 

z  z    2   +        + 

To  evaluate  the  effect  of  this  operator  on  the  wave  functions  (61) ,  it  is 

necessary  to  recall  that  the  action  of  the  operators  L  ,  L  and  L_  on  a 

wave  function  |  £ni  ^>  give 
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(65) 


Lz  |  £m£>    =  m£  |  ta^, 

L+l£rV>    =Q(*+1)   -  m£(m£+l^|jlmA+l> 

Similarly,    for   later  reference,    the   spin  operators   S    ,    S      and  S_,    acting 
on  a   function    |sm  ^S  ,    result  in 

S    Ism  ^>=  m    Ism  "^ 

z '      s-^         s '      s<^ 

(66) 

S    |Sm  ^>=rs(S+l)    -  m    (m  +l7|     Ism  +1^>    . 
+ '       s  *^       L_  s      s—  _J      '       s 


Operating  with  Lz  on    (61)    gives  Ljl^Sm  ^>  =   /2    (|  33-m/>  +    |3-3-m\V 

The  matrix  element  of  L     using  this  basis   is   then 

z 

<r2Sms|Ljr2Sms>=  |<33|ms|33|ms>-  §  <3-3|ma  |  3-3|mg>=  0. 

The  matrix  elements  of  L+  and  L_  using  the  Equations  (61)  as  basis  functions 
are  also  zero,  since  the  function  (61)  does  not  contain  terms  which  have 
mi   ±1*   The^efore  the  matrix  elements  vanish  by  orthogonality.   The  first 
term  of  (62)  is  therefore  zero. 

The  third  term  of  (62)  is  evaluated  by  using  the  identity 

L-T  =  L  S   +  -j  (L+S_  +  L_S+)  .  (67) 

The  effect  of  operating  on  | r  S- ^  with  L-s\  for  example,  is  then 

r-n  r2s|>  =  ^[  3X|  i>  -J  33f  |> 

"i  1 3-4  !> + 1 3-4  !>  •  <m> 
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Operating  from  the  left  on  Equation  (68)  with  each  of  the  four  wave  func- 

3    1 
tions  (61)  formed  by  using  mg  =  ±-,  ±-  gives  zero  in  every  case.   In  a 

similar  way,  it  is  easy  to  show  that  the  matrix  elements  of  "iMs  are  all 

zero  for  the  other  basis  functions  (61).   The  third  term  of  Equation  (62) 

is  therefore  zero,  showing  that  the  fourfold  spin  degeneracy  of  the  ground 

state  is  not  lifted  by  the  first  order  spin-orbit  coupling  perturbation. 

The  only  non-zero  matrix  elements  of  the  first  order  perturbation  are  then 

<Cr2Smslgs6H-s|r2Sms^>  .  (69) 

The  energy  contribution  from  the  second  order  perturbation  is 


^      |  <^  T2Sms  I  0-  (t  +   gt)    +  hl-t\  TSms  ^> 


(70) 


r  =  r  r  E(r)  "  E(roJ 

V  5 


where  ?2   represents  the  ground  state  and  r  represents  the  excited  states 

r4'F5*   The  ener9ies  of  these  states  are  E(T  )  and  E(D,  respectively. 
The  term  i 

<T„Sm  |g  "sl  TSm  "^> 
2   s '  s  '    s^ 

of  the  expression  (70)  is  zero  because  the  operator  S*  affects  only  the 

Sms  parts  of  the  wave  functions,  leaving  the  orbital  term  <^r  1 1^>  which 

vanishes  because  T      and  T   are  orthogonal.   The  second  order  perturbation 


then  becomes 


I 


<^r2sm  |bh-T+  AL.1T|rsm  ^>|2 

S •  (71) 


P  E(T)    -   E(r    ) 


4'. 


The  numerator  of  the  above  expression  is  expanded  to  give 
|\r2Sms|AL-"^|rsms^|   +  2  <^T   Sm  |  B&-"l"|  FSm  ^> 

•  ^rsm  lAL-s|r  Sm  *S  +  I  <J\Sm  |etl-"fc|rsm  ^>  I2.     (72) 
^   s      '2   s*"^      ^*  2   s         s^ 

The  last  term  of  (72)  does  not  depend  on  spin,  therefore  it  shifts  each  of 

the  four  ground  state  levels  by  the  same  amount.   Since  this  term  causes  no 

26 
change  in  relative  ground  state  energy,  it  is  usually  ignored.    Combining 

Equations  (69) ,  (70)  and  (72)  gives  the  ground  state  energies  correct  to 

second  order  of  the  perturbation  (60)  as 

Er2  =  <r2Smsl^*lr2Sms> 


I 


,  E(T)   -  E(r2) 


<r  Sm  |AL-"s|rsm  ^> 
2      s1  '  s-^ 


+  2  <r2sms|6H.L|rsms>  (73) 

•  <TrSm  |AL-"s|r  Sm  ^> 
^    s1     '2   s-^J 

As  noted  earlier,  the  first  term  of  (73)  is  the  only  contribution  to  the 

first  order  energy  of  the  T      states.   If  all  other  terms  were  zero,  this 

term  alone  would  cause  an  isotropic  splitting,  linear  with  magnetic  field, 

of  the  fourfold  spin  degenerate  T      level  into  four  separate  equally  spaced 

sublevels  characterized  by  m  =  ±-  and  m  =  ±-.   The  second  term  in  the 

1      s     2       s     2 

summation  accounts  for  the  anisotropy  of  the  PMR  spectrum,  and  shows  that 
this  anistropy  is  due  to  second  order  spin-orbit  coupling  between  the  V 
and  T  states.   The  first  term  in  the  summation  gives  the  splitting  of  the 
ground  state  in  the  absence  of  a  magnetic  field.   When  a  magnetic  field  is 
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applied,  each  of  the  two  doublet  levels  is  further  split  into  two  sub- 
levels.   Transitions  between  these  various  levels  are  called  the  PMR  fine 
structure.   A  graph  of  this  splitting  for  YAG:Cr  is  included  in  Section  V. 

Another  contribution  to  the  energy  (73)  which  has  been  ignored  is 
the  contribution  due  to  the  nuclear  magnetic  moment  of  the  active  ion. 
If  one  of  the  isotopes  of  this  ion  has  half-odd  nuclear  spin  I,  the  pos- 
sible spin  states  of  the  nucleus  will  cause  21+1  incremental  values  of 
magnetic  field  to  be  superimposed  upon  the  external  magnetic  field,  thus 
splitting  the  PMR  absorption  for  this  isotope  into  21+1  separate  ab- 
sorptions.  This  phenomenon  is  called  the  nuclear  hyperfine  interaction, 
and  its  effect  on  the  PMR  absorption  is  called  the  hyperfine  structure. 
Chromium  has  one  isotope  (53)  with  nuclear  spin  |  and  9.5  percent  abun- 
dance, and  one  isotope  (54)  with  integral  nuclear  spin.   Because  of  the 
low  percentage  of  chromium  53  in  ordinary  samples,  the  hyperfine  struc- 
ture is  difficult  to  observe  and  has  not  been  observed  in  this  work.   For 

this  reason,  the  nuclear  spin  contributions  to  the  energy  (73)  have  been 

27 
ignored.   Bleaney  and  Bowers   have  observed  chromium  hyperfine  structure 

in  dilute  chrome  alums  enriched  with  chromium  53. 

III.B.   Determination  of  the  Ground  State  Energy  Levels 
^—  — —u. 

with  H  Parallel  to  the  Symmetry  Axis 

If  the  external  magnetic  field  is  applied  parallel  to  the  crystal- 
line symmetry  axis,  the  four  ground  state  levels  will  diverge  linearly, 
and  the  g-f actor ,  which  is  defined  for  this  orientation,  may  be  found  by 
measuring  the  splittings  of  the  m  =  ±-  levels.   The  zero  field  splitting 
is  also  found  in  this  orientation  by  extrapolating  the  splittings  measured 
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at  high  fields  to  zero  field.   If  the  magnetic  field  is  applied  at  an 
angle  to  the  symmetry  axis ,  there  is  competition  between  the  crystal  field 
and  the  magnetic  field  to  determine  the  axis  of  precession  of  the  spinning 
electrons,   and  the  energy  levels  do  not  diverge  linearly.   This  competition 
results  from  the  fact  that  the  ground  state  has  a  small  angular  momentum 
due  to  the  combined  action  of  the  trigonal  distortion  and  the  second  order 
spin-orbit  interaction. 

Values  of  both  the  g-factor  and  the  zero  field  splitting  are  required 
to  complete  the  characterization  of  the  crystalline  field  begun  in  the  last 
section.   For  this  reason,  Equation  (73)  will  be  used  to  calculate  the 
energy  levels  with  the  magnetic  field  and  symmetry  axes  aligned,  because 
this  approach  provides  insight  into  the  way  the  ground  state  behaves  when 

perturbed.   Other  orientations  will  be  treated  by  using  the  spin  Hamiltonian 

28  29 

formalism  devised  by  Pryce   and  Abragam  and  Pryce  "  as  applied  to  chromium 

9  .  . 

in  octahedral  symmetry  by  Davis  and  Strandberg.    This  treatment  was  origi- 
nally applied  to  ammonium  chrome  alum,  but  it  also  accurately  predicts  the 
locations  of  PMR  transitions  in  YAG:Cr. 

The  basis  functions  which  will  be  used  for  the  calculation  of  ground 

state  energies  are  those  given  by  Equations  (33)  and  (33a)  for  E  =  -12Dq 

and  E  =  -2Dq  and  (42) ,  (43)  for  E  =  2Dq  and  E  =  lODq.   The  result  of  this 

calculation  will  then  be  used  to  determine  the  r    level  splitting  which 

0  ~~2 
will  yield  the  trigonal  distortion  parameter  A  r  .   By  including  the 

0  2 
effects  of  A  r  ,  a  new  set  of  crystal  field  eigenvalues  and  eigenfunctions 

can  be  calculated,  as  was  done  in  Section  II,  the  results  of  which  are  given 
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in  Table  I.   These  new  eigenfunctions  will  then  be  used  to  re-evaluate 
the  ground  state  splittings  to  include  any  effects  of  the  trigonal  dis- 
tortion.  In  particular,  it  will  be  of  interest  to  evaluate  the  effects 

4 
of  the  admixed   P  states . 

If  the  magnetic  field  is  aligned  with  the  crystalline  symmetry  axis, 
Equation  (73)  becomes 

Er   =  <J\Sm  k   H  s  lr-,Sm  ^ 
T_    ^  2   s   s   z  z'  2   s 


T   e(d  -\(T2)  I    I  <r2sms|At.t|rsms>  I 


2 


♦  2<r2Smsl6H2Lz|rsms>  (74) 

•  <JSmjALZS2|r2Sms^ 


•>] 


1  3 

The  matrix  elements  of  g  BH  S   are  simply  ±-g  BH  and  ±-g  BH  .   The  operator 

S  Z       Z  Z       S  Z  A       S  Z 

BH  L     couples  only  states    | 1^    and    | 11^   and   its  matrix  element   is   indepen- 
dent of  spin.    Operating  on   the   state    1 1^   with  L     gives 

L    |I>    =^L    |30>-4l    |33>    +^L    |3-3> 
z'«^  3      z '     -^  3      z'     ■y  3      z1       -y 

=  -»/2|33>    -    /2|3-3> 

Operating  with   \£l|    gives 

<ll|Lz|l^>   =   2,  (75) 

and  all  of  the  other^T    \l,    |  f^  are  zero.   Operating  on  the  states 

1 1  —  +  — ^  with  L  S  gives 
'   2  —  2"^       z  z 
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v.iifif>-4v.i«l±£>-4v>4±£> 


3      z    z '         2  -  2^  3      z   z '         2  -  2' 

4lS    I  3-3  f  -  |> 

3      z    z  2   —  2^ 


+  —  L   S 


J33|±|>--||3-3|±|> 


II  ,  +  fl  gives 
<ZI  f  ±  flLzSz  I1  1±1>=+3-  <76> 


In  a  similar  way,  the  other  matrix  elements  of  L  S   are  found  to  be 

z  z 

^ItflVj1!!^-^-  (76a, 

The  matrix  elements  of  L-S  are  found  by  operating  on  each  of  the  states 
lr2Sms^with  L*S  as  was  done  in  arriving  at  Equation  (76).   The  result- 
ing functions  are  then  operated  on  with  all  of  the  bra  vectors*^  Sm  I 

^  4   s  ' 

and  \r5msl-   These  matrix  elements  are  listed  in  Table  II. 


TABLE  II 


2  s- 

-3 


Matrix  Elements  ^Tm  Il-sIt  m  ^ 
>  s1      2  s^ 

<«l|M|if> 
<mi|M|i|>.  ft 

<ni|t.?|ii>.    -! 
<III-i|t.?li>=  zfi 

<iv||t-S|ii>   =,  JS 

<H-i|t.«|l-i>    .   ! 


(continued) 
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TABLE   II    (continued) 

<m-||M|i-i>-  G 
<ivi|M|i-i>  =  2S2 
<ii-§|M|i-j>   =  3 

<lV-i|M|l-§>    =  /6 

Note  that  all  matrix  elements  between  r  and  T     states  are  zero.   This 

•  4   4 

result  is  also  obtained  if  the  F/  P  interaction  is  not  considered  in 

deriving  the  eigenvalues  and  eigenvectors  of  the  crystal  field. 

The  results  of  Equations  (75)  ,  (76)  and  Table  II  are  now  substituted 
into  Equation  (74)  to  obtain  the  ground  state  energies.   These  energies 


are : 


IcBH-A'TK^lMliOl2 


Ero3  =  2  gsBH  "  A       '^    2' Li 

22  I  En "  EI 


|<mi|M|i|> 


E  -    E 

III  I 


<Il||LS     |l|> 
\     21    z    z1    2^ 


2A3H 

EirEi 


OlEJ"> 


(77a) 


3  2 

=  I  g  BH  -  a" 
2      s 


12A3H 


r  9         6   1 

[En  -  ei  +  Em  -  eiJ 


En  -  ei  ' 


r       =  i  g  BH  -  a2   I"— +  " +  — -1 

Y2\        2      S  LEH    "   EI        Em    "   EI        EIV  "   ElJ 


E 
1", 

4A3H 


E         -   E    ' 
II  I 


(77b) 
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E       x  =  -\  gs3H  -  A2 
2 


r^ — + — § — + ~j 1 

LEn _  Ei  Em  ■  Ei  eiv  -  eiJ 


4AgH 

+  f ZT~  '  <77c) 

EII    EI 


E    3  =  -\    gs3H  -  A2 
2 


["-i s_1 

I  E    -  E     E     -  E 
|_  II     I     III     ij 


L  12AgH 

+  s — rj-  •  (7?d) 

ii      i 

The  zero-field  splitting  is  obtained  by  setting  H  =  0  in  Equations  (77)  , 
and  results  in 

Er      ,   "  Er      n    =   6   =  8A2 


F2+|  r2+i 

-2  -2 


r—i i— i 

E  -   E  E         -   E 

Lin       i       ii       ij 


where  the  fact  that  E.^  =  E   is  used.  This  energy  difference  is  given 

to  good  approximation  by 

8A2(En  -  Em) 
6  =  — -  2 '  (79) 


(\  "  \) 


where  Er  is  the  average  value  of  the  energy  of  the  rc  level  and  E  = 

r5  5  Z 

E   .   The  effective  g- factor  is  obtained  by  subtracting  (77c)  from  (77b) 
2 

to  get 

8A3H 


2       2  5     2 

Using  the  condition  for  observation  of  PMR  gives 

AE  =  hv 


•91  -«(«.-»). 


5     2 
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so  that 

8A 

9  =  9_  -  p ~S —  *  (80) 

s    E     E 

5      2 

This  equation  shows  that  the  free  electron  g-factor  is  modified  by 

second  order  spin-orbit  interaction  between  T„  and  T_  states.   It  will 

z       b 

be  shown  in  Section  V  that  Equation  (80) ,  together  with  the  measured 

values  of  g  and  E   -  E   ,  give  a  value  of  A  that  yields  poor  results 

5     2 

when  used  in  the  subsequent  crystal  field  calculation.   For  this  reason, 
the  free  ion  value  91  cm   quoted  by  Low   will  be  used  for  determi- 
nation of  the  trigonal  distortion  coefficient. 

If  the  proper  value  of  A  is  used,  Equation  (79)  yields  the  T     level 

splitting  in  terms  of  the  ground  state  splitting  6,  determined  from  PMR 

measurements,  and  the  position  of  the  r_  level  Er  -  Er  ,  determined 

5        F5    F2 

from  optical  absorption  measurements. 

III.C.   Determination  of  the  Ground  State  Energy 
Levels  for  an  Arbitrary  Magnetic  Field  Orientation 

The  spin  Hamiltonian  is  widely  used  to  determine  the  ground  state 
energy  levels  of  paramagnetic  ions.   This  approach  reduces  the  orbital  de- 
pendance  of  the  perturbation  (60)  to  constant  numbers  that  have  only  the 
function  of  parameters,  and  it  is  then  only  necessary  to  evaluate  the  matrix 
elements  of  the  spin  operators  to  obtain  the  spin-Hamiltonian  matrix  ele- 
ments.  The  spin  Hamiltonian  therefore  operates  only  on  spin  states,  and 

its  matrix  elements  and  energy  levels  are  evaluated  by  using  these  states. 

28  . 
Pryce   has  expressed  the  spin  Hamiltonian  in  the  elegant  form 
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Hs  =  e(g^6^  "  2AAi^  S'H-  "  a2a-  -s-s-  (81) 

in  which  the  summation  convention  is  employed.   In  Equation  (81)  ,  s  repre- 
sents spin,  H  is  magnetic  field,  and  the  quantity  A.,  is  the  tensor 

A     _?<rAl£><rhlr2> 

ij  •      E(T)  -  E  (T  ) 


Note  the  similarity  of  Equation  (81)  to  the  analogous  equation  (73). 

29 
Abragam  and  Pryce   have  expanded  the  spin  Hamiltonian  formalism  to  in- 
clude the  effects  of  the  nuclear  hyperfine  interaction,  which  need  not 
be  considered  for  YAG:Cr. 

The  spin  Hamiltonian  approach  has  been  applied  to  chromium  in  octa- 
hedral symmetry  by  Davis  and  Strandberg,9  who  used  it  to  explain  PMR  in 
ammonium  chrome  alum.   They  have  solved  the  4  x  4  secular  determinant, 
corresponding  to  mg  =  +-,  +-,  in  terms  of  a  normalized  energy  and  mag- 
netic field  W  =  2E/6  and  X  =  2g3H/6,  where  6  is  the  ground  state  split- 
ting.  Their  solution  to  the  secular  determinant  is 

4     2       5   2  I       ?  ? 

W   -W    2  +  2XJ   +  2X  W  Pcos  8  -  1)  +  1 

_„2     2n        1   2     9   4 
-  3X  cos  0  +  -  X  +  —  X  =0,  (82) 

where  9  is  the  angle  which  the  magnetic  field  makes  with  the  crystalline 
symmetry  axis.   This  equation  is  solved  in  Section  V,  using  measured  values 
of  g  and  6,  for  several  different  orientations,  and  the  energy  levels  will 
be  plotted  as  a  function  of  magnetic  field. 
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III.D.   Calculation  of  Ground  State  Energy  Levels 
Including  Trigonal  Distortion  Effects 

The  calculation  of  energy  levels  carried  out  in  Section  III.B.  was 
repeated  using  the  wave  functions  of  Table  I,  which  include  the  trigonal 
distortion  effects.   To  perforin  this  calculation,  the  corrected  wave 
functions  were  substituted  into  Equation  (74)  to  give  expressions  analo- 
gous to  Equations  (77)  for  the  ground  state  energies. 

As  before,  the  matrix  elements  of  g  3H  S  are  simply  ±-   a   BH  and 

szz  2   s   z 

3 
2  gs^Hz"   The  matrix  elements  of  L  are  obtained  by  first  operating  on 

the  wave  function  1 1 ' '*>  with  L  to  obtain 

"^  z 

^l1' C>  =   "(3)  (0.457)  |33>   -(3)  (0.457)  |  3-3>    . 

The  only  wave   function  which  couples   to  the  above  expression  is   the   func- 
tion   | II' '^>  ,    so  that 

<rlLzlr2>=    <II"ILZII">   =   "    1-94.  (83) 

The  matrix  elements  of  L  S     are  given  by 


<5Ii't*s«lIi>-  - 2-91 


<J4K*J4>--°- 


968 


(84) 


The  other  matrix  elements  of  L-T  are  determined  in  the  same  way  as  the 
matrix  elements  given  in  Table  II,  but  these  new  elements  will  not  be 
given  here.   Using  the  above  results,  the  parameter  6  is  given  to  good 
approximation  by 

7.54A2(ElI  -  ElII) 
6  *  (E         -  E  )2  '  (85) 
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and  the  corrected  g- factor  is 

7.78A 
5     l2 
These  results  show  that  the  trigonal  distortion  has  little  effect  on  the 
parameters  6  and  g,  as  might  be  expected  in  view  of  the  slight  changes  in 
state  function  coefficients  caused  by  the  distortion.   Of  more  significance 
is  the  fact  that  all  matrix  elements  of  L-S  between  T      and  T      states  no 
longer  vanish.   Most  of  these  elements  remain  small  enough  to  have  negli- 
gible effects  on  the  energy  levels,  but  the  L-S  interaction  will  still  mix 
small  amounts  of  higher  lying  states  with  the  T      ground  state.   This  mixing 
probably  accounts  for  the  observation  of  PMR  transitions  which  do  not  obey 

the  magnetic  dipole  transition  selection  rule,  both  in  this  work  and  in  the 

9      •»  •*• 
work  of  Davis  and  Strandberg.    The  L-S  perturbation  simply  mixes  a  higher 

lying  state  of  the  proper  spin  with  the  ground  state  in  such  a  way  that  the 

selection  rule  is  satisfied  for  a  small  component  state  of  the  wave  function 

and  a  weak  transition  is  therefore  observed. 


SECTION  IV 
EXPERIMENTAL  PROCEDURE 

IV. A.   Sample  Preparation 

The  sample  of  YAG:Cr  studied  during  this  work  was  obtained  from  the 
Airtron  Division  of  Litton  Industries.   The  sample  was  finished  in  the 
form  of  a  6  mm  x  12  mm  cylinder  and  was  specified  to  have  0.1  percent  of 
the  aluminum  ions  replaced  by  chromium  ions.   At  this  dopant  level,  the 
crystal  is  very  pale  green  —  almost  colorless.   This  green  color  is 
caused  by  broad  optical  absorption  bands  centered  at  wavelengths  of  431 
and  595  nanometers.   Although  the  dopant  level  was  specified  to  be  0.1 
percent,  no  attempt  was  made  to  verify  this  concentration. 

The  sample  was  mounted  in  a  goniometer,  which  was  then  placed  in  a 
back  reflection  Laue  X-ray  diffraction  camera.   A  Laue  photograph  was 
made  of  the  crystal  which  showed  that  the  cylindrical  axis  was  very 
nearly  parallel  to  the  LlllJ  direction.   The  goniometer  was  adjusted  to 
precisely  align  the  x-ray  beam  with  the  [ill}  direction,  and  another 
photograph  was  made  to  verify  this  alignment.   The  goniometer  was  then 
removed  from  the  camera,  with  this  alignment  preserved,  and  a  face  was 
ground  on  the  crystal,  parallel  to  the  (111)  planes,  with  a  diamond  abra- 
sive wheel.   This  procedure  was  repeated  for  the  QLlcQ  and  [poll    direc- 
tions, resulting  in  a  crystal  with  faces  ground  parallel  to  the  (111), 
(110)   and  (001)  sets  of  planes.   Figures  3a,  3b,  and  3c  show  Laue  photo- 
graphs made  along  each  of  these  directions.   Note  that  the  spot  pattern 
on  each  photograph  displays  the  characteristic  symmetry  of  the  axis  along 
which  the  crystal  was  oriented  for  the  photograph.   Each  of  these  pic- 
tures was  made  using  a  molybdenum  x-ray  tube  operated  at  19.5  kilovolts 
and  30  milliamperes ,  and  was  exposed  for  30  minutes. 
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Figure  3a.  Laue  X-ray  Diffraction  Pattern 
of  YAG:Cr  Sample  in  the  QlllJ  Direction 
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Figure  3b.  Laue   x-ray  Diffraction  Pattern 
of  YAG:Cr  Sample  in  the  Qioj  Direction 
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Figure  3c.   Laue  X-ray  Diffraction  Pattern 
of  YAG:Cr  Sample  in  the  foo^  Direction 
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The  end  of  the  crystal  on  which  the  faces  were  ground  was  cut  off 
and  retained  for  PMR  measurements.   The  remaining  piece,  a  cylinder  of 
dimensions  about  6  mm  x  8  mm,  was  polished  on  both  ends  for  optical  ab- 
sorption measurements. 

IV. B.   Optical  Absorption  Measurements 

The  optical  absorption  spectrum  of  the  YAG:Cr  sample  was  measured 
on  a  Beckman  Model  DK-2  spectrophotometer.   This  instrument  is  used  in 
many  laboratories  to  measure  optical  transmission  of  filters,  laser 
mirrors,  chemicals,  and  other  optical  materials  in  applications  where 
extreme  accuracy  is  not  required.   It  is  suitable  for  measurements  on 
YAG:Cr  because  the  absorption  bands  are  broad  and  ill-defined,  so  that 
use  of  a  more  sophisticated  instrument  is  not  required. 

The  Beckman  DK-2  is  a  dual  channel  prism  spectrometer  which  uses 
the  second  channel  for  reference  in  cancelling  the  effects  on  the  spectrum 
due  to  source  and  detector  variations  with  wavelength.   It  uses  a  tung- 
sten-iodide quartz  lamp  as  a  source,  and  a  photomultiplier  tube  or  lead 
sulfide  cell  as  a  detector,  depending  on  the  wavelengths  to  be  studied. 
The  DK-2  covers  a  wavelength  range  of  200  to  4000  nanometers.   The  re- 
sults of  optical  absorption  measurements  made  on  YAG:Cr,  using  the  DK-2, 
will  be  discussed  in  Section  V. 

IV. C.   Paramagnetic  Resonance  Measurements 

Microwave  spectrometers  for  the  detection  of  PMR  absorptions  were 
designed  and  built  in  the  course  of  this  work  for  measurements  at  fre- 
quencies near  10   Hz  (X-band)  and  3.5  x  10   Hz  (Ka-band).   These 

30  31  32 
spectrometers  are  of  conventional  design,  '   '   except  that  they  do  not 
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use  the  superheterodyne  detection  technique  used  in  the  most  sensitive 
spectrometers.  The  sensitivity  of  these  spectrometers  is  adequate  for 
the  observation  of  the  relatively  strong  absorptions  which  occur  in  YAG:Cr. 

Figure  4  is  a  block  diagram  of  the  10  GHz  system.   Microwave  energy,, 
obtained  from  an  X-band  klystron,  is  coupled  to  the  system  through  an 
isolator  and  E-H  tuner.   A  small  amount  of  the  energy  is  picked  off  at 
the  cross  guide  coupler  and  mixed  with  a  harmonic  from  a  1  MHz  frequency 
standard  and  fed  back  to  lock  the  klystron  to  the  standard  for  stable 
operation.   The  radio  frequency  oscillator  mixes  with  the  frequency 
multiplier  output  to  provide  frequencies  between  the  harmonics  of  the 
standard.   The  majority  of  the  klystron  energy  is  coupled  through  an 
attenuator  to  one  port  of  a  magic  tee.   This  device  divides  the  energy 
equally  between  two  outputs,  one  of  which  goes  to  the  tuned  microwave 
cavity  between  the  poles  of  an  electromagnet,  and  the  other  goes  to  an 
amplitude  and  phase  altering  system  consisting  of  an  attenuator  and  a 
slotted  guide  tuner.   This  system  is  adjusted  until  the  energy  reflected 
from  it  just  balances  the  energy  reflected  from  the  cavity  in  which  the 
sample  is  mounted.   Therefore,  little  or  no  energy  leaves  the  fourth 
port  of  the  magic  tee  unless  an  absorption  occurs  in  the  cavity.   This 
bridge  arrangement  minimizes  the  quiescent  signal  incident  on  the  detec- 
tor and,  therefore,  reduces  detector  noise.   The  sample  cavity  is  tuned 
to  resonate  at  the  applied  microwave  frequency,  with  a  cavity  mode  con- 
figuration chosen  such  that  the  microwave  and  external  magnetic  fields 
are  orthogonal.   This  resonant  condition  is  characterized  by  a  high 
cavity  "Q"  or  quality  factor,  and  a  sharply  defined  absorption  of  micro- 
wave energy  at  the  resonant  frequency.   When  a  PMR  absorption  occurs  in 
the  sample,  microwave  energy  is  absorbed  and  the  Q  of  the  cavity  changes, 
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thus  increasing  the  amount  of  signal  incident  on  the  fourth  port  of  the 
magic  tee.   Energy  from  this  fourth  port  is  then  detected  and  fed  into  a 
lock-in  amplifier,  which  drives  a  strip  chart  recorder  on  which  the  PMR 
spectrum  is  recorded.   The  reference  output  of  the  lock-in  is  used  to 
drive  modulation  coils  mounted  on  the  pole  pieces  of  the  magnet.   Fre- 
quency measurement  is  accomplished  by  picking  off  part  of  the  energy  from 
the  adjustable  arm  of  the  magic  tee,  zero-beating  it  with  a  harmonic 
of  the  transfer  oscillator,  and  counting  the  transfer  oscillator  fre- 
quency.  Precise  magnetic  field  measurements  are  made  with  a  nuclear 
magnetic  resonance  gaussmeter.   The  cavity  is  suspended  in  a  double 
dewar  and  may  be  cooled  if  necessary,  but  cooling  was  not  required  to 
observe  the  strong  resonances  in  YAG:Cr.   In  this  system,  the  microwave 
frequency  is  held  constant  while  the  magnetic  field  is  swept  slowly. 
PMR  absorptions  appear  as  differentiated  lines  because  of  the  method  of 
modulation  and  detection  used. 

Figure  5  is  a  block  diagram  of  the  35  GHz  PMR  spectrometer.   The 
operation  of  this  system  is  almost  identical  to  the  operation  of  the 
10  GHz  system,  and  will  not  be  discussed  in  detail  except  for  the  few 
important  differences  between  the  two  systems.   The  10  GHz  phase  locked 
source  discussed  earlier  is  reduced  in  frequency  to  less  than  9  GHz  and 
is  used  as  a  stable  reference  for  the  35  GHz  spectrometer,  which  is  phase 
locked  to  the  fourth  harmonic  of  this  source.   The  output  frequency  of 
the  35  GHz  system  is  then  a  factor  of  four  higher  than  the  frequency  of 
the  9  GHz  input,  plus  or  minus  the  10  MHz  reference  frequency,  and  can 
be  precisely  determined  by  measuring  the  frequency  of  the  9  GHz  source 
with  the  transfer  oscillator  and  frequency  counter.   A  10  MHz  ambiguity 
in  frequency,  introduced  by  the  reference,  can  be  easily  resolved  by 
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using  the  wave  meter  to  give  a  coarse  frequency  measurement.   A  duplexer 
is  used  in  place  of  the  magic  tee  for  coupling  energy  to  the  sample 
cavity.   This  device  is  a  ferrite  coupler  used  for  directing  transmitted 
and  received  energy  to  and  from  antennas  in  microwave  radars  and  com- 
munications systems. 

The  sample  cavities  for  both  spectrometers  were  fabricated  by  using 
an  electroforming  process  in  which  a  stainless  steel  mandrel  of  the  de- 
sired dimensions  was  machined  and  polished  to  a  very  smooth  finish.  A 
thin  layer  of  gold  was  then  deposited  on  the  mandrel,  followed  by  the 
main  cavity  deposition  which  was  about  2  mm  of  copper.   The  mandrel  was 
then  pulled  from  the  electro formed  deposition,  leaving  a  cavity  with 
inner  surfaces  gold  plated  and  with  very  precisely  defined  corners. 
Flanges  were  then  soldered  to  the  cavity  and  the  proper  microwave  coup- 
ling iris  sizes  were  determined  by  increasing  the  hole  size  until  the 
desired  degree  of  coupling  was  obtained. 

The  cavity  used  for  10  GHz  measurements  was  rectangular  and  the  35 
GHz  cavity  was  circular  in  cross  section.   Because  of  space  limitations 
in  the  dewar,  the  10  GHz  cavity  was  fix-tuned  and  would  support  only  one 
resonant  mode.   The  35  GHz  cavity  was  tunable  and  was  large  enough,  be- 
cause of  the  higher  frequency,  to  support  several  modes.   Tuning  was  ac- 
complished by  a  movable  plunger  coupled  by  stainless  steel  rods  to  a 
micrometer  screw  mounted  on  the  dewar  flange. 

For  observation  of  PMR  absorptions,  the  external  magnetic  field  and 
the  microwave  magnetic  field  must  be  orthogonal.   In  a  single  mode  rect- 
angular cavity,  this  orthogonality  may  be  achieved  simply  by  mounting 
the  sample  centered  in  the  bottom  of  the  cavity.   However,  if  the  external 
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magnetic  field  is  rotated  to  obtain  data  on  other  orientations,  this 
orthogonality  is  lost  with  a  resultant  reduction  in  signal  strength. 
This  problem  is  solved  by  mounting  the  sample  on  the  side  wall  of  the 
cavity  as  shown  in  Figure  6a,  with  the  result  that  the  two  fields  are 
orthogonal  regardless  of  the  orientation  of  the  electromagnet.   The 
sample  may  be  mounted  centered  on  the  bottom  of  the  cylindrical  cavity, 
because  a  cavity  mode  can  be  found  in  which  the  magnetic  field  is  axial, 
as  shown  in  Figure  6b. 

Figure  7  shows  how  PMR  in  several  crystal  orientations  of  interest 
can  be  observed  in  a  cubic  crystal  without  the  necessity  for  removing 
and  remounting  the  sample.   Faces  are  ground  on  the  crystal  parallel  to 

(111),  (110),  and  (001)  planes  as  described  in  Section  IV. A.   The  cry- 
stal is  then  mounted  in  the  cavity  with  its  (110)  planes  horizontal,  so 
that  the  external  field  is  always  parallel  to  these  planes  when  the 
field  is  rotated.   The  external  field  may  then  be  aligned  with  the  [ill] , 

[lioj  ,  and  j_00lj  directions,  as  well  as  many  others,  as  the  external 
electromagnet  is  rotated,  without  remounting  the  sample. 


Sample , 


Sample 
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Figure  6.   Methods  of  Mounting  Samples  in  Rectangular  and  Circular  Cavities 

so  That  the  External  and  Microwave  Magnetic  Fields  Are  Always  Orthogonal. 

The  External  Field  May  Be  Oriented  in  Any  Direction  in  a  Horizontal 

Plane  Perpendicular  to  the  Plane  of  the  Page. 
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Figure  7.   Orientation  and  Rotation  of  Sample  for  PMR  Measurements. 


SECTION  V 
RESULTS 

V.A.   Paramagnetic  Resonance  Measurements 
PMR  absorption  spectra  were  measured  at  both  X-band  and  Ka-band. 
The  X-band  measurements  were  made  as  the  magnetic  field  was  rotated  in 
10   steps  relative  to  the  crystalline  symmetry  axis,  beginning  with  the 
field  parallel  to  the  Lliqj  direction  and  ending  with  the  field  parallel 
to  the  JOOlj  direction,  as  shown  in  Figure  7.   The  Ka-band  measurements  were 
made  with  the  magnetic  field  parallel  to  the  [lioj,  [ill],  and  Q)0l] 
directions  only,  because  these  directions  are  considered  most  important. 
Actually,  only  the  measurements  made  in  the  H-ll]  direction  are  neces- 
sary for  characterization  of  the  trigonal  distortion  of  the  crystal 
field. 

Figures  8,  9,  and  10  show  the  results  of  Ka-band  measurements  made 
for  the  three  orientations  of  interest.   The  X-band  measurements  are  not 
shown  because  they  are  not  as  easily  interpreted  as  the  Ka-band  results. 
These  lower  frequency  results  are  complicated  by  the  fact  that  the 
energy  of  the  zero-field  splitting  is  greater  than  the  energy  of  an  X- 
band  photon,  resulting  in  a  large  grouping  of  transitions  at  lower  fields 
where  lines  are  weaker  and  may  be  broadened  by  ill-defined  energy  levels. 

In  a  cubic  crystal,  there  are  four  possible  equivalent  0-11]  direc- 
tions.  If  an  external  magnetic  field  is  aligned  with  one  of  these  direc- 
tions, the  field  direction  will  make  an  angle  0  with  that  direction  and 
an  angle  70.6  with  the  other  three  directions.   Measurement  of  PMR  ab- 
sorptions along  this  particular  direction  should  then  show  contributions 
from  the  Cr   sites  displaced  parallel  to  the  field  as  well  as  those 
displaced  at  angles  of  70.6   to  the  field  direction.   This  assertion  is 
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Figure  8.   Paramagnetic  Resonance  in  YAG:Cr  [ill]  Direction,  v  =  35.73  x  10S 
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Figure  9.  Paramagnetic  Resonance  Spectrum  of  YAG:Cr  [_lloJ  Direction,  v  =  35.73  x  10  Hz 
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Figure  10.   Paramagnetic  Resonance 


Spectrum  of  YAG:Cr  [ool]  Direction,  v  =  35.73  GHz 
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confirmed  by  the  PMR  spectrum  for  the  Ql^  direction  shown  in  Figure  8. 
Note  that  the  lines  are  generally  grouped  into  triplets  and  singlets. 
The  singlet  lines  are  due  to  ions  displaced  parallel  to  the  field,  and 
the  triplets  are  contributed  by  ions  displaced  at  an  angle  of  70.6   to 
the  field.   Ideally,  the  triplets  should  be  degenerate,  but  it  is  pos- 
sible that  some  unexplained  crystal  field  anomaly,  or  a  slight  misorien- 
tation  of  the  crystal,  has  caused  them  to  be  split.   The  center  line  of 

the  triplet  located  at  12.895  kilogauss  is  actually  degenerate  with  the 

o  11 

0  line  which  results  from  the  transition  m  =  +-  •*""*—-. 

s    2     2 

If  the  external  magnetic  field  is  aligned  parallel  to  the  crystal- 
line |_-'--'-0j  direction,  two  of  the  trigonal  displacement  directions  make 
an  angle  of  90  with  the  field  and  the  remaining  two  make  an  angle  of 
35.3  with  the  field.   One  would  then  expect  the  lines  to  occur  in  doub- 
lets, an  assertion  which  is  confirmed  in  some  measure  by  Figure  9.  Note 
that  some  of  the  doublets  are  degenerate  while  others  are  split  by  more 
than  100  gauss.   When  the  external  field  is  aligned  with  the  LpOlJ  direc- 
tion, it  makes  an  angle  of  54.7  with  each  of  the  trigonal  displacements, 
so  that  the  absorption  lines  should  occur  in  quartets,  which  is  confirmed 
by  Figure  10. 

The  zero  field  splitting  is  calculated  by  using  the  PMR  transitions 
measured  in  the  [ill]  direction  which  are  attributed  to  ions  displaced 
parallel  to  the  magnetic  field.   A  more  accurate  calculation  of  this 
parameter  is  possible  if  transitions  are  used  which  have  high  field 
values,  because  the  locations  of  these  transitions  are  known  to  better 
accuracy.   This  higher  accuracy  is  a  result  of  the  splittings  varying 
almost  linearly  with  magnetic  field  at  the  higher  field  levels.   To  make 
this  calculation,  the  locations  of  two  Ka-band  and  one  X-band  transition 
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were  used,  as  shown  in  the  energy  level  schematic  of  Figure  11 


Figure  11.   Ground  State  Energy  Level  Schematic 
with  External  Magnetic  Field 


By  using  this  schematic,  it  is  not  difficult  to  show  that  the  zero  field 
splitting  6  is 


6  = 


H2\  -  Vx 
H1+H2 


(87) 


where  Ek  and  Ex  are  energies  of  Ka  and  X-band  photons,  respectively, 
and  the  magnetic  field  intensities  H  ,  H  are  defined  by  the  figure. 
The  desired  Ka-band  resonance  occurs  at  a  frequency  of  35.73  GHz  and  a 
magnetic  field  of  7.377  kilogauss,  while  the  X-band  resonance  occurs  at 
10.071  GHz  and  9.354  kilogauss.   Substituting  these  results  into  (87) 
gives  6  =  0.509  +  0.001  cm   ,  where  the  basis  for  the  accuracy  will  be 
discussed  later.   The  effective  g-factor  is  calculated  from 


hVk  =  Ek  =  g3H' 


(88) 
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where  the  variables  have  been  defined  earlier.   For  this  calculation,  it 

is  necessary  to  use  the  parameters  corresponding  to  the  transition 

1     1  9 

ms  =  ~2  **-  +2>    namely  H  =  12.895  kilogauss,  and  v  =  35.73  x  10  Hz.   This 

calculation  gives  g  =  1.980  ±  0.002.   These  results  agree  well  with  those 

of  Carson  and  White   who  found  6  =  0.510  cm"1  and  g  =  1.98,  for  YAG  doped 

with  chromium. 

Using  the  above  values  for  6  and  g,  it  is  now  possible  to  solve 
Equation  (82)  for  the  normalized  energy  level  splittings  as  a  function  of 
angle  and  magnetic  field.   This  equation  was  solved  for  the  [ill]  ,  HlIoI  , 
and  [poi]  directions  and  the  results  are  plotted  in  Figures  12,  13,  and  14. 
The  PMR  transitions  which  have  been  observed  in  the  course  of  this  experiment, 
both  at  X-band  and  at  Ka-band,  are  designated  by  arrows.   The  long  arrows 
indicate  Ka-band  transitions  and  the  short  arrows  represent  X-band  tran- 
sitions.  Equation  (82)  was  used  to  compare  the  locations  of  measured  tran- 
sitions observed  in  the  [ill]  direction  to  the  locations  predicted  by  the 
theory  of  Davis  and  Strandberg.   The  observed  values  of  magnetic  field  were 
substituted  into  (82)  and  this  equation  was  solved  for  the  energies,  which 
were  then  compared  to  the  microwave  energy  measured  for  the  transitions. 
The  results  of  these  calculations  are  given  in  Table  III. 

TABLE  III 
PMR  Line  Positions  in  the  Q.ll]  Direction  at  35.73  GHz 


Cr   Ion 


Magnetic  Field  in  Gauss 


Angular 

Displacement  Calculated  Measured  Percent  Difference 

0°  3,445  3,533  2.55 

0°  7,371  7,377  0.08 

0°  12,895  12,895  0.00 

70.6°  5,238  5,239  0.02 

(continued) 
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TABLE  III  (continued) 


Cr  +  Ion 

Magnetic  Field 

in  Gauss 

Angular 

Displacement 

Calculated 

Measured 

Percent  Difference 

70.6° 

6,954 

6,933 

-0.30 

70.6° 

10,637 

10,691 

0.51 

70.6° 

12,913 

12,895 

-0.14 

70.6° 

14,373 

14,369 

-0.03 

V.B.   Optical  Absorption  Measurements 

The  optical  absorption  spectrum  of  YAG:Cr  was  measured  over  the 
spectral  range  250  to  900  nanometers,  using  the  Beckman  DK-2  spectro- 
photometer, and  the  results  are  shown  in  Figure  15.  No  additional  ab- 
sorptions were  observed  when  the  measurement  range  was  extended  from 
200  to  2800  nanometers.   This  spectrum  shows  three  broad  absorptions, 
whose  positions  and  qualitative  widths  are  predicted  with  good  accuracy 
by  crystal  field  theory.   The  accuracy  of  the  measured  energies  shown 
in  the  figure  is  estimated  to  be  +  100  cm  .   Although  it  is  possible 
to  achieve  much  higher  accuracy  by  using  a  more  sophisticated  spectro- 
meter, the  uncertainty  in  absorption  position  caused  by  the  widths  of  the 
absorptions  does  not  justify  this  additional  accuracy. 

The  absorptions  shown  in  Figure  15  are  broadened  by  the  trigonal 
distortion  splitting,  spin-orbit  coupling  within  the  manifold,  and  by 
thermal  vibration  of  the  ions  which  give  rise  to  the  crystal  field.   Un- 
fortunately, the  latter  two  effects  effectively  mask  the  trigonal  dis- 
tortion splitting,  so  that  this  effect  could  not  be  measured  by  optical 
absorption  techniques,  and  the  value  of  the  splitting  calculated  in  the 
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Figure  12.   Normalized  Ground  State  Splitting  Versus 
Magnetic  Field  for  Qll]  Direction 
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Figure  13.   Normalized  Ground  State  Splitting  Versus 
Magnetic  Field  for  QllOJ  Direction 
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Magnetic  Field   in  Kilogauss 

Figure    14.      Normalized  Ground  State  Splitting  Versus 
Magnetic  Field   for    jTlOl]    Direction 
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Figure  15.   Optical  Absorption  Spectrum  of  YAG  Doped  with  0.1%  Chromium 
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next  section  cannot  be  verified  experimentally.   Some  degree  of  experimental 
verification  is  possible  by  comparing  the  calculated  splitting  to  the 
measured  half -width  of  the  optical  absorption,  a  comparison  which  is  made 
in  Figure  17. 

V.C.   Determination  of  the  Crystal  Field  Parameters 

■ ■  - 

In  Section  II,  it  was  shown  that  the  crystal  field  potential,  including 
the  trigonal  distortion,  may  be  written  in  the  form 


r  °  <  ,v/2 , 3  -37i 


V=A°r2y°+D-    r4 


,02         4 
where  A2  r   and  D'  r  are  parameters  to  be  determined  by  PMR  and  optical 

absorption  measurements.   It  was  proven  in  Section  III  that  the  splitting 

of  the  T5  crystal  field  level,  in  terms  of  parameters  defined  earlier,  is 

&(e 


(%  -  \) 2 


EII  "  Eiu  ■  A  "  ' 2 "*"  '  (90> 

8A 

This  equation  gives  A  in  terms  of  measurable  quantities  if  the  spin-orbit 

coupling  constant  A  is  known.   It  is  tempting  to  use  Equation  (80) ,  which 

gives  g  in  terms  of  A  and  E   -  E    to  calculate  A  based  on  PMR  measure- 

5     2 
ments.   Unfortunately,  this  procedure  gives  a  value  for  A,  which  when  sub- 
stituted into  (90)  ,  results  in  a  value  of  A  that  is  not  physically  reason- 
able.  Using  the  values  gg  =  2.0023,  g  =  1.980  from  PMR  measurements,  and 
A  =  16,700  cm  '  from  optical  absorption  measurements,  Equation  (80)  gives 
A  =  46  cm   .   Substitution  of  this  value  of  A  into  (90)  results  in  A  = 
8400  cm   which  is  more  than  one-half  of  the  separation  of  the  T     and  r 
states.   This  value  of  A  does  not  agree  with  the  half-width  of  the  T 
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energy  level  determined  by  optical  absorption  measurements  and  is  therefore 
not  physically  reasonable.   This  lack  of  agreement  can  probably  be  blamed 
on  covalency  effects,  studied  by  Owen,    and  the  need  for  correcting  the 
free  electron  g-f actor  for  the  atomic  environment,  discussed  in  detail  by 


33 
Abragam  and  Van  Vleck. 


The  problem  of  choosing  the  proper  value  of  A  for  this  work  has  been 
solved  by  using  the  free  ion  value.   Kleiner   gives  theoretical  and  experi- 
mental evidence  that  the  value  of  the  spin-orbit  coupling  parameter  A  has 
a  negligible  dependence  on  the  crystalline  field.   He  shows  that  A  may  have 
a  contribution  AQ,  valid  for  the  free  ion,  and  a  contribution  A  due  to  the 


crystalline  field.   He  then  shows  that  A  /A  %   10~4  so  that  the  free  ion 
value  AQ  is  a  good  approximation  to  A  in  a  crystalline  field.   The  fact 
that  magnetic  susceptibility  data  for  several  iron  group  ions  agree  well 
with  theory  based  on  the  free  ion  spin-orbit  coupling  parameter  also  supports 

this  contention.   The  appropriate  free  ion  value,  consistent  with  treating 

■*••*■  10-1 

L-S  as  a  state  operator,  is  given  by  Low,    and  is  91.0  cm   .   Substituting 

this  value  of  A  and  the  previously  derived  value  of  6  into  the  equation  for 
A  gives  A  =  (2.14  ±  0.03)  x  103  cm"1. 

It  must  be  noted  that  the  sign  of  A,  the  T      level  splitting,  is 

ambiguous.   Equation  (79)  shows  that  if  A  =  E   -  E    is  positive,  then 

the  ground  level  splitting  6  is  positive  in  the  sense  that  E   +3  is  greater 

F2  "2 
than  E   +1  so  that  A  and  6  have  the  same  sign.   By  using  optical  and  PMR 

2  ~2 
absorption  measurements,  it  is  unfortunately  not  possible  to  determine  the 

sign  of  either  6   or  A.   It  is ,  however,  possible  to  show  that  A  must  be 
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positive  by  solving  the  cubic  equation  (55)  for  A°  r2  using  both  positive 
and  negative  values  of  A,  and  thereby  showing  that  a  positive  value  for  A 
makes  the  most  physical  sense.   Substituting  the  negative  value  for  A, 
together  with  the  value  Dq  =  1670  cm"1,  into  (55)  gives  A°  r2  =  -  2.61  x 
10   cm   ,  while  the  positive  value  of  A  gives  1.33  x  104  cm"1,  3.38  x  104  cm"1, 
and  2.79  x  10   cm   .   In  the  former  case,  the  other  two  roots  of  the  cubic 
equation  are  complex,  and  therefore  not  physically  meaningful.   The  proper 

value  of  A2  r   is  chosen  by  substituting  A  =  2140  cm"1  into  the  equation 

for  the  T5  level  splitting  correct  to  first  order,  and  obtaining  A°  r2  = 

15,900  cm   .   The  root  of  Equation  (55)  nearest  this  value  is  1.33  x  104 

obtained  by  substituting  the  positive  value  of  A  into  (55).   Considering 

the  indicated  tolerances  for  Dq  and  A,  it  is  therefore  concluded  that 

0  ~2  4         -1 

A  r  =  1.33  x  10   ±  200  cm   . 

The  crystal  field  parameter  Dq  is  determined  from  the  location  of 
the  T  level  as  given  by  optical  absorption  measurements.   Figure  15 
shows  the  optical  absorption  spectrum  of  YAG:Cr.   The  energy  of  the  V 
level  is  simply  lODq  as  shown  in  Section  II,  so  that  Dq  =  1670  +  10  cm 
It  is  reasonable  to  ask  whether  the  second  order  trigonal  distortion 
perturbation  will  affect  the  T     and  I"  levels  enough  to  make  it  neces- 
sary to  correct  Dq  for  this  effect.   Equations  (49) ,  (54)  and  (56)  are 

used  to  evaluate  these  corrections.   Substituting  the  above  value  of 

0  ~2  .         . 
A  r   into  Equation  (56)  gives  for  the  corrected  T      energy 

E   =  -12Dq  -  160  cm"  . 


5 
-1 
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Similarly,  Equation  (56)  gives  for  the  energy  E   of  the  T      level 

E    =  -2Dq  +  1200  cm"1, 
and  Equation  54  gives  for  the  E    and  E   levels 

EII3.  =  E    =  -2Dq  -  950  cm"1. 
Considering  the  twofold  degeneracy  of  the  E    and  E   states,  the  aver- 
age energy  of  the  r _  level  is  then 

D 

Ep   =  -2Dq  -  235  can"1. 

5 

These  results  show  that  the  energy  difference  E   -  E_  is  lODq  -  75  cm~ 

5     2 
which  is  essentially  lODq  within  experimental  error.   Based  on  these 

calculations,  it  is  concluded  that  the  second  order  trigonal  distortion 

correction  has  negligible  effect  on  the  determination  of  Dq. 

— 4 
The  coefficient  D'r  of  the  fourth  order  term  in  the  crystal  field 

expansion  may  now  be  calculated  using  the  above  determined  value  of  Dq 

and  one  of  the  diagonal  crystal  field  matrix  elements  determined  in 

Section  II.   Using  the  3,3  matrix  element,  this  calculation  gives 

|  <33|y°|33>  D'r4  =  -2Dq, 


or 


1  (^Aiiij  D'r4  =  ~2^ 


which  gives 

D'r4  =  -  4.14  x  104  +  250  cm"1. 

2 
The  factor  -  in  the  above  equation  results  from  expressing  the  crystal 

field  potential  in  a  form  such  that  the  [ill]  direction  is  the  symmetry 
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axis,  as  noted  in  Section  II.   Substituting  for  A°  r2  and  D'  r4  in  the 
crystal  field  perturbation  expansion  gives  finally 

V  =  13,300  y°  -  41,100  [y°  +  (1°,  ^  (yJ  -  y^H   ,       (91) 

expressed  in  cm"  . 

V.D.   Accuracy  of  Results 

Errors  in  frequency  measurement,  magnetic  field  measurement,  and 

crystal  orientation  are  the  types  of  errors  which  must  be  considered 

when  evaluating  the  accuracy  of  PMR  measurements.   A  study  of  these 

3     1 
errors  was  made  in  measuring  the  parameters  of  the  -  •*— ►  -  transition 

in  YAG:Cr  at  X-band.   Both  frequency  and  magnetic  field  for  this  tran- 
sition were  measured  a  total  of  10  times  while  keeping  the  magnetic 
field  centered  on  the  absorption  line.   After  each  measurement,  the 
field  was  readjusted  to  the  center  of  the  line  and  the  frequency  and 
magnetic  field  were  measured  again.   During  this  series  of  measurements, 
the  frequency  remained  stable  at  10.7067  GHz,  while  the  magnetic  field 
varied  from  9353.7  to  9354.6  gauss.   The  mean  of  the  magnetic  field 
measurements  was  9354.1  gauss  and  the  standard  deviation  was  0.28  gauss. 
Assuming  that  a  random  field  measurement  will  fall  within  3  standard 
deviations  (3a)  of  the  mean,  the  magnetic  field  for  this  transition  is 
9354.1  +0.8  gauss  or  the  accuracy  is  +  0.009  percent.   This  accuracy 

percentage  is  assumed  for  all  of  the  magnetic  field  measurements,  and 
the  frequency  measurement  errors  are  considered  negligible,  based  on 
the  measured  stability  of  the  frequency  source. 


8(> 


Perhaps  the  largest  error  contribution  in  the  PMR  measurements  is 
caused  by  a  slight  misorientation  of  the  crystal  relative  to  the  mag- 
netic field  direction.   This  misorientation  may  be  caused  by  errors  in 
X-ray  orientation,  crystal  mounting,  or  spectrometer  alignment.   The 
Laue  patterns  shown  in  Figures  3a,  3b,  and  3c  show  that  the  X-ray  orien- 
tation is  probably  accurate  to  a  few  tenths  of  a  degree  because  a  1 
error  in  direction  will  cause  about  a  2  mm  deflection  of  the  central  X- 
ray  spot.   The  photographs  show  that  the  central  spots  are  centered  in 
the  pattern  within  less  than  about  0.5  mm.   The  standard  deviation  of 
alignment  error  caused  by  X-ray  orientation  is  then  estimated  to  be  0.3 
degrees.   The  crystal  was  oriented  in  the  PMR  spectrometer  by  aligning 
the  intersections  of  the  faces  ground  on  the  crystal  with  the  microwave 
cavity  geometry  using  a  surveyor's  telescope.   The  standard  deviation 
of  this  error  is  estimated  to  be  0.5  degrees.   Finally,  the  spectrometer 
was  aligned  with  the  electromagnet  axis  by  again  using  the  surveyor's 
telescope.   Because  of  the  large  sizes  of  the  magnet  and  spectrometer, 
this  error  is  considered  negligible  compared  to  the  other  two  error 
sources.   Assuming  that  the  alignment  errors  are  normally  distributed, 


[2        2l  l 
(0.5)   +  (0.3)H 

=0.58  degrees.   One  would  then  expect  the  alignment  error  to  be  almost 

always  less  than  (3) (0.58)  =  1.74  degrees. 

To  evaluate  the  effect  of  the  above  error  on  the  location  of  the 
PMR  resonances,  Equation  (82)  was  solved  using  the  parameters  determined 
for  the  m  =  --  *—*■   -  transition  at  35  GHz  for  values  of  8  ranging  from 
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-5°  to  +5°.   The  results  of  this  calculation  are  plotted  in  Figure  16. 
Note  that  an  alignment  error  of  1.7  degrees  corresponds  to  an  error  in 
magnetic  field  determination  of  13  gauss,  which  is  large  compared  to  the 
field  measurement  error  and  is  therefore  the  dominant  error  source.   This 
magnetic  field  error  was  substituted  into  Equations  (87)  and  (88)  for  6 
and  g  to  give  the  error  limits  indicated  in  the  evaluation  of  these 
parameters.   The  error  in  6,  together  with  the  error  in  optical  absorp- 
tion measurements,  was  then  used  to  estimate  the  uncertainty  in  the  Tr 
level  splitting,  A. 

The  3a  error  for  the  optical  absorption  measurements  is  estimated 
to  be  100  cm   .   Error  contributions  in  these  measurements  come  from 
spectrometer  calibration  alignment  as  well  as  uncertainty  in  position  of 
the  absorption  maximum  due  to  the  width  of  the  absorption  itself.   This 
error  in  the  location  of  the  optical  absorption  maxima  was  used  to  esti- 
mate the  error  in  determination  of  the  crystal  field  parameters  Dq  and 
D'  r  . 

V.E.   Summary  and  Conclusions 
The  results  of  this  study  are  summarized  in  Table  IV.   In  those  in- 
stances where  previously  known  results  are  available,  they  are  given  for 
comparison.   Figure  17  is  a  series  of  energy  level  diagrams  based  on  the 
calculations  and  measurements  of  the  parameters  given  in  Table  IV.   This 
Figure  shows  the  effect  of  the  first  and  second  order  trigonal  distor- 
tion perturbation  on  the  energy  levels,  as  well  as  the  results  obtained 
by  solving  the  complete  10  x  10  crystal  field  matrix  by  computer.   The 
first  energy  level  diagram  shown  is  that  measured  by  optical  absorption, 
and  the  shaded  regions  are  intended  to  represent  the  absorption  half- 
widths.   Absorption  peaks  are  indicated  in  these  shaded  regions.   The 
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percentage  errors  between  measured  and  calculated  centers  of  gravity  are 
indicated  on  the  figure.   Degeneracies  are  noted  in  parentheses. 

TABLE  IV 
Summary  of  Results 

Value 
This  Work         Carson  and  White 


Parameter 

Crystal  Field  Parameter, 
Dq  (Measured) 

Crystal  Field  Expansion 
(Calculated) 


Ground  State  Splitting, 
6  (Measured) 

r5  State  Splitting,  A 
(Calculated) 

Spectroscopic  Splitting 
Factor  g  (Measured) 


Location  of  Ta   Levels 
(Measured) 


Location  of  I"  Levels 
(Calculated) 


1670   +  10  cm 


-1 


0  |~o 

2-41,400[y4 

3 


-i 


13,300y 

h. 
,,10.     ,    3      -3 

+(— )   (y4-y4 

0.509  +  0.001  c 
2.14  x  103 
1.980  +  0.002 


2.32  x  104  +  100  cm  1 
3.57  x  104  +  100  cm"1 


2.37  x  104  cm-1 
3.71  x  104  cm-1 


0.510 


1.98 


The  most  important  result  of  this  study  is  the  use  of  complementary 
optical  and  paramagnetic  resonance  spectroscopy  to  characterize  the  crys- 
talline field  in  YAG:Cr.   The  link  between  these  two  different  spectro- 
scopic methods  is  the  connection  between  the  trigonal  distortion  split- 
ting of  optical  levels  and  the  slight  splitting  of  the  ground  level  which 
has  such  a  large  effect  on  the  PMR  spectrum. 
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Figure  17.   Calculated  and  Measured  Energy  Level  Diagrams 
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The  locations  of  the  optical  absorption  bands  due  to  4F  and  P 
multiplet  states  have  been  successfully  predicted  by  using  crystal  field 
theory.   The  widths  of  these  bands  have  been  qualitatively  accounted  for 
by  treating  the  trigonal  distortion  as  a  perturbation  of  the  octahedral 
crystal  field  energy  levels.   The  locations  of  PMR  absorptions  have  been 
found  to  agree  well  with  the  theory  of  Davis  and  Strandberg. 

The  method  of  treating  the  trigonal  distortion  as  a  perturbation  of 
the  octahedral  crystal  field  has  application  to  other  ions  in  other  cry- 
stalline complexes.   For  example,  the  neodymium  ion  in  YAG  occupies  a 
site  of  orthorhombic  but  near  cubic  symmetry.   It  is  possible  that  the 
deviation  from  cubic  symmetry  could  be  treated  as  a  perturbation  of  the 
cubic  states  to  better  explain  the  energy  levels  of  this  important  laser 
material.   In  addition,  there  are  many  materials  containing  iron  group 
ions  to  which  this  procedure  might  be  applied  to  better  explain  observed 
magnetic  resonance  and  optical  absorption  results. 
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